NASA  TM  X-  72800 


NASA  TECHNICAL 

MEMORANDUM 


NASA  TM  X-  72800 


AD-A279  566 


FLUTTER  OF  LAMINATED  PLATES 
IN  SUPERSONIC  FLOW 
BY 


JAMES  WAYNE  SAWYER 


This  informal  documentation  medium  is  used  to  provide  accel  erated  or 
special  release  of  technical  information  to  selected  users.  The  contents 
may  not  meet  NASA  formal  editing  and  publication  standards,  may  be  re¬ 
vised,  or  may  be  incorporated  in  another  publication. 


NATIONAL  AERONAUTICS  AND  SPACE  ADMINISTRATION 
UNGLEY  RESEARCH  aNTER,  HAMPTON,  VIRGINIA  23665 


94- ‘14434 

IIIUIIin'TIIIMIIIIIIIII 


1.  Report  No.  2.  Gowimmtnt  Acconion  No. 

NASA  TM  X- 72800 

3.  Recipient's  CMaiog  No. 

4.  Title  and  Subtitle 

Flutter  of  Laminated  Plates  in  Supersonic  Flow 

S.  Raport  Data 

Nov.  1975 

6.  Performing  Organizatiop  Coda 

7.  Authorial 

James  Wayne  Sawyer 

8.  Performing  Organization  Raport  No. 

10.  Work  Unit  No. 

506-17-22-01 

9.  Informing  Organitation  Nama  and  Addrau 

NASA  Langley  Research  Center 

Hampton,  VA  23665 

11.  Contract  or  Grant  No. 

13.  Type  of  Report  and  Pariod  Covered 

Technical  Memorandum 

12.  Sponsoring  Agency  Name  and  Address 

National  Aeronautics  and  Space  Administration 
Washington,  DC  20546 

14.  Sponsoring  Agency  Code 

16.  Supplementary  Notes 

16.  Abstract 

A  solution  procedure  has  been  developed  using  linear  small  deflection  theory 
for  the  flutter  of  simply  supported  laminated  plates.  For  such  plates,  the 

bending  and  extenslonal  governing  equations  are  courled  and  have  cross¬ 
stiffness  terms  which  do  not  appear  in  classical  place  theory.  An  extended 
Galerkin  method  is  used  to  obtain  approximate  solutions  to  the  governing 
equations,  and  the  aerodynamic  pressure  loading  used  in  the  analysis  is  that 
given  by  linear  piston  theory  with  flow  at  arbitrary  cross-flow  angle.  A 
limited  parametric  study  has  been  conducted  for  typical  laminated  composite 
plates.  The  calculations  show  that  both  the  bending-extensional  coupling  and 
the  cross-stiffness  terms  have  a  large  destabilizing  effect  on  flutter. 

Since  classical  plate  theory  does  not  consider  bending-extensional  coupling 
and  cross  stiffness  terms,  it  usually  gives  inaccurate  and  nonconservative 
flutter  boundaries  for  laminated  plates. 


17.  Kay  Words  (Suggastad  by  AuthorlsM  (STAR  category  underlined! 

18.  Distribution  Statamant 

STRUCTURAL  MECHANICS 

Panel  Flutter 

Composite  Plates 
Laminated  Plates 

Unclassified 

-  Unlimited 

19.  Sacurity  Oasaif.  (of  this  report) 

30.  Sacurity  ClasMf.  (of  this  page) 

31.  No.  of  PagM 

22.  Wee' 

Unclassified 

Unclassified 

119 

.i;5.25 

Thi  National  Technical  Information  Sirvice,  Springfield,  Virginia  32151 

'Available  from 

STIF/NASA  Scientific  and  Technical  Information  Facility,  P.O.  Box  33,  Collaga  Park.  MO  30740 


TABLE  OF  CONTENTS 


CHAPTER  PAGE 

TITLE .  ^ 

TA3LE  OF  CONTENTS .  11 


1. 


2. 


NOTATION .  Iv 

LIST  OF  TABLES .  vii 

LIST  OF  FIGURES .  viii 

INTRODUCTION  .  1 

1.1  Review  of  Literature  .  2 

1.1.1  Formulation  of  Equations  .  3 


1.1.2  Solutions  for  Symmetric  Plates  . 

1.1.3  Solutions  for  Angle-ply  and  Cross-ply  Plates  . 


1.1.1*  Approximate  Solutions .  8 

1.1.5  Dynamic  Stability  .  9 

1.2  Statement  of  Problem  .  9 

DEVELOPMENT  OF  GOVERNING  EQUATIONS  .  13 


2.1  Lamina  Stress-strain  Equations  .... 

2.2  Lamina  and  Plate  Constitutive  Equations 

2.3  Plate  Governing  Equations  . 

2.3.1  Equations  of  Motion  .  .  . 

2.3.2  Compatibility  Equation  . 

2.h  Boundary  Conditions  . 


ii 


.  .  .  1*^ 

.  .  .  16 

21 


21 


26 

28 


By_ 

Distribution^ 
Availability  6o4mi 

retail  ai»5/or~' 

Spaolal 


/  i 

fci  or  : 

•id 

'.loa 


OO 


lii 

CHAPTER  PAGE 

3.  APPROXIMATE  SOLUTIONS  OF  THE  FLUTTER  EQUATIONS  .  31 

3.1  General  Laminated  Plates  .  32 

3.2  Symmetric  Laminated  Plates .  ^1 

3.3  Angle-ply  Plates . 

3.^  Approximate  Solutions . ^9 

U.  RESULTS  AND  DISCUSSION . 52 

l*.l  Convergence  of  Results .  53 

h,2  Comparison  with  Literature . 5^ 

k.3  Symmetric  Plates . 56 

U.3.1  Flutter  Boundaries  .  5T 

U.3.2  Effects  of  Inplane  Normal  and  Shear  Loads.  .  56 

U.3.3  Effect  of  Cross- flow  . .  59 

U.U  Angle-ply  Plates  .  6l 

4.1*.l  Flutter  Boundaries . 6l 

U.U.2  Effect  of  Inplane  Normal  and  Shear  Loads  .  .  63 

U.U.3  Effect  of  Cross-flow . 63 

U.5  General  Laminated  Plates  .....  . 

U.3.1  Flutter  Boxindaries . 63 

U.5.2  Effect  of  Plate  Construction  .  65 

U.5.3  Composite  Stiffened  Aliuninum  Plates  ....  6T 

5.  CONCLUSIONS  AND  RECOMMENDATIONS . 69 

REFERENCES . 


NOTATION 


C  ,  C 
mn*  rs 


^11*^22 


H  ,H 
inn  rs 


K  ,  K 
x’  y 


matrix  of  elastic  areas  defined  by  eqviatlon  (2.l8) 

matrix  defined  by  equation  (2.23) 

coefficient  matrix  (See  equation  (3.18).) 

plate  dimensions  along  x  and  y  eocls,  respectively 

matrix  of  elastic  statical  moment  defined  by  equation  (2.l8) 

matrix  defined  by  equation  (2.23) 

coefficient  matrix  (See  equation  (3.18).) 

Fotirier  series  coefficients 
coefficient  matrix  (See  eqxiation  3.18.) 

matrix  of  elastic  moments  of  inertia  defined  by  equation  (2.18) 

matrix  defined  by  equation  (2.23) 

matrix  of  with  6*0. 

matrix  of  with  9“0  and  6'=0 

Principal  Young's  moduli  in  1  and  2  directions,  respectively 
Airy  stress  function,  defined  by  equation  (2.31) 
shear  modulus  in  1-2  planes 
totcLl  plate  thickness 

dlsteoice  from  plate  midplane  to  the  lamina  boundary  as  shown 
in  figure  3. 

Fourier  series  coefficients 
integers 


number  of  layers  in  plate 
nondimen.sion  inplane  normal  loads 


iv 


V 


K 


xy 


00 

M,N 

M  ,M 
X  y 

M 

xy 

m,n 


N  ,  N 
x’  y 


nondimensional  inplane  shear  load. 


2 

a 


2 

IT 


Mach  number 

toted  number  of  terms  in  x  emd  y  direction,  respectively 
bending  moment  resultants  as  defined  by  equation  (2.9) 
twisting  moment  resulteint  as  defined  by  equation  (2.9) 
integers 

inplane  normal  loads  per  unit  length  in  x  and  y 
directions ,  respectively 


N' 


xy 


N 

xy 

N'  ,  N' 

X  y 

N  ,  N 
x’  y 


N 


xy 


«ij 


QiJ 


<1 


induced  inplane  shearing  loads 
inplane  shearing  loads  per  unit  length 
induced  inplane  normal  loads 

inplane  normal  resultant  as  defined  by  equation  (2.8) 
inplane  shear  resvdtant  as  defined  by  equation  (2.8) 
intensity  of  transverse  load 

constitutive  coefficient  matrix  for  specially 
orthotropic  lamina 

constitutive  coefficient  matrix  for  generally  orthotropic 
lamina 

dynamic  pressure  of  airstream 
integers 

midplane  displacements  in  x  and  y  direction,  respectively 
normal  displacements 


4 


x,y,z 


cartesian  coordinates 


vl 


6  compressibility  factor* 

Y  plate  mass  per  unit  area 

Y  >Y  >Y  shearing  strain  components 

xz 

5  lamina  thickness  (See  figure  1^) 


normal  strains  in  x,  and  y  directions*  respectively 
nomal  midplane  stredns  in  x  and  y  directions*  respectively 
rotation  angle  of  the  fibers  with  respect  to  the  plate  axis 


9'  rotation  angle  of  the  fibers  in  the  inner  laminas  with  respect 

to  plate  axis  (See  figiires  ll»  euid  15.)  • 

*^x’'^y*'^xy  plate  curvature  as  defined  by  equation  (2.5) 

A  cross-flow  angle 

X  flutter  parameter  defined  by  equation  (3.16) 

poisson's  ratio  for  orthotropic  lamina  defined  by  -  ®2^®i 

0*0*0  normal  stress  components  in  x*  y,  and  z  directions,  respectively 
X  y  z 

T  *  T  ,  T„_  shearing  stress  components 
xy  xz’  yz  ^ 

(1)  natural  frequency 


u 


r 


fundamental  frequency  of  simply  supported  beam  as  given 
by  equation  (3.17) 


LIST  OF  TABLES 


TABLE  PAGE 

I.  LAMINA  MATERIAL  PROPERTIES  USED  IN  ANALYSIS . 78 

II.  MATERIAL  PROPERTIES  FOR  COMPOSITE  STIFFENED  ALUMINUM  PLATE.  .  79 

III.  MATERIAL  PROPERTIES  OF  REFERENCE  PLATE  .  80 

IV.  COMPARISON  OF  FOUR  LOWEST  NATURAL  FREQUENCIES  FOR  A  SQUARE 

ANGLE-PLY  PLATE  ^12  *  '  .  .  .  8l 

V.  COMPARISON  OF  FOUR  LOWEST  NATURAL  FREQUENCIES  FOR  A 

SQUARE  CROSS-PLY  PLATE  (G^g/E^^  =  .5,  =  .25) . 82 


vli 


LIST  OF  FIGURES 


FIGURE  PAGE 

1.  -  Geometry  of  angle-ply  plates . 83 

2.  -  Geometry  of  cross-ply  plates . 81* 

3.  -  Lamina  coordinate  system . 85 

U.  -  Lamina  notation  and  positive  direction  of  stress  .  86 

5<  -  Rotation  of  lamina  axes .  8? 

6.  -  Positive  direction  of  stress  and  moment  resviltants .  88 

7.  -  Plate  geometry  suid  coordinate  system .  89 

8.  -  Nomenclature  and  geometry  of  symmetric  plates .  90 

9.  -  Nomenclature  and  geometry  of  angle-ply  plates . 91 

10.  -  Nomenclature  and  geometry  of  general  laminated  plates  ....  92 

11.  -  Nomenclatiire  and  geometry  of  composite  stiffened  aluminum 

plate .  93 

12.  -  Convergence  of  solution  for  a  symmetric,  boron-epoxy,  square 

plate  with  K  =  U .  93 

13.  -  CosuLescence  of  frequencies  for  a  symmetric,  boron-epoxy, 

squaire  plate  with  K=U .  96 

lU.  -  Comparison  of  flutter  solutions  for  square  reference  plate. 

N  =  M=  2 .  98 


15. 


16. 


17. 

18. 


Converged  flutter  boundaries  for  square  reference  plate  .  .  . 

Flutter  boundaries  for  square,  symmetric  plates.  N  =  N  = 

N  =0 . *..... 

xy 

Flutter  boundaries  forjboron-epoxy,  symmetric  plate  with 

a/b  =  2.0.  N  =N  =N  =0  . 

X  y  xy 

Flutter  boundaries  for  square,  boron-epoxy,  symmetric  plate 
with  inplane  normal  loads.  N  =N  =0;K=U  . 

y 


99 

100 

102 

103 


vill 


lx 


FIGURE  PAGE 


19.  -  Flutter  boundaries  for  sqviare,  boron-epoxy,  symmetric  plate 

with  inplane  shear  loads,  if  sjT  =  0;K=U . lOU 

y  X 

20.  -  Coalescence  of  frequencies  for  a  symmetric,  boron-epoxy, 

square  plate  with  inplane  shear  loads.  N  =N  =0;K=U; 

and  e  =  15° . f . 105 

21.  -  Flutter  boundaries  for  square,  boron-epoxy,  symmetric  plate 

with  cross-flow,  n  =n  =N  =0;K  =  U . IO8 


X  y  xy  ’ 

22.  -  Flutter  boundaries  for  square,  angle-ply  plate. 

N  =  N  =  N  =0 . 109 

X  y  xy 

23.  -  Flutter  boundaries  for  boron-epoxy,  angle-ply  plate  with 

a/b  =  2.0.  N  =  N  =  N  =0 . Ill 

x  y  xy 

2U.  -  Flutter  boundaries  for  squeire,  boron-epoxy,  angle-ply  plate 

with  inplane  loads.  K  =  U . 112 


25.  -  Flutter  boundaries  for  square,  boron-epoxy,  angle-ply  plate 

with  cross-flow.  N  =N  =  N  s0;K  =  l4 . llU 

X  y  xy  ’ 

26.  -  Flutter  boundaries  for  square,  boron-epoxy,  general 

laminated  plates.  N=N=N^  =0 . 115 

X  y  xy 

27.  -  Comparison  of  flutter  boxmdaries  for  square,  symmetric, 

angle-ply,  and  general  laminated  plates,  N  =  N  =  N  =  0  .  117 

X  y  xy 

28.  -  Effects  of  cross-flow  on  flutter  of  square,  symmetric, 

angle-ply  and  general  laminated  plates.  N  =  N  =  N  =  0.  .  II8 

X  y  JQT 

29.  -  Flutter  boundaries  for  square,  composite-stiffened  aluminvm 

plates.  N  =N  =N  =0;  mass  ,  .  /mass  ..  .  ^ 

^  X  y  xy  aluminum  composite  =  1.0 


119 


Chapter  1 


INTRODUCTION 


Laminated  composite  materials  have  been  in  use  for  many  centuries. 
Pieces  of  laminated  wood  have  been  found  which  date  back  to  about  1500  B.C. 
Some  excellent  swords  were  made  in  the  15th  century  by  laminating  severad 
layers  of  steel  to  provide  an  extremely  hard  and  keen  cutting  edge  with 
a  softer  tough  body.  Although  some  of  the  advantages  of  using  composite 
materials  have  been  known  for  some  time,  only  in  recent  years  have  deter¬ 
mined  efforts  been  made  to  fully  develop  composite  materials.  The  advent 
of  high  performance  aircraft  and  spacecraft  has  brought  about  the  need 
for  more  efficient  use  of  materials.  Recent  energy  shortages  have  em¬ 
phasized  in  even  a  greater  way  the  importance  of  developing  composite 
materials  for  everyday  usage. 

In  recent  years,  composite  materials  have  been  the  subject  of  a  large 
number  of  experimental  emd  analytical  investigations  and  have  been  develop¬ 
ed  to  the  point  where  structural  components  are  now  being  used  in  actual 
applications.  Previously,  most  of  the  composite  applications  have  been 
structured  components  built  to  replace  existing  conventional  components. 
However,  new  designs  are  emerging  which  have  considered  composite  materials 
in  the  total  design.  As  new  applications  are  considered  for  composite 
designs,  a  better  understanding  of  the  behavior  of  structural  components 
under  various  load  conditions  Is  needed. 

One  structviral  component  that  is  used  extensively  in  a  wide  variety 
of  applications  is  the  flat  plate.  Plates  fabricated  from  composite 
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materials  usually  consist  of  individual  lamina  bonded  together.  Each 
lamina  usuedly  consists  of  a  nxjmber  of  high  strength  filaments  aligned  in 
the  same  direction  aind  held  in  place  by  a  plastic  matrix  material. 

The  lamina  may  be  highly  orthotropic ,  and  since  the  principal  material  direc¬ 
tion  or  the  angle  the  fibers  make  with  the  plate  axis  may  differ  for  each 
lamina,  the  resulting  anisotropic  plate  is  more  complex  to  analyze  than 
an  orthotropic  plate. 

This  study  is  being  conducted  to  obtain  a  better  understanding  of 
the  flutter  characteristics  of  laminated  composite  flat  plates.  However, 
the  following  discussion  will  include  a  general  review  of  previous  studies 
of  laminated  composite  plates.  Most  previous  studies  of  laminated  plates 
have  been  concerned  only  with  deflection,  vibration,  and  buckling  char¬ 
acteristics.  The  development  and  solution  of  the  governing  equations 
for  plates  under  various  conditions  will  be  reviewed  and  a  proposed  area 
of  investigation  identified.  Solution  techniques  and  possible  problem 
areas  will  also  be  discussed. 

1.1  Review  of  Literature 

Nxmerous  theoretical  analyses  of  composite  laminated  plates  have  been 
made  and  are  discussed  in  references  1-23.  These  references  represent 
work  that  cover  the  spectrum  from  classical  smadl  deflection  theory  and 
approximate  solutions  to  non-linear  theory  and  exact  solutions.  Much  of 
the  early  analytical  work  for  composite  plates  was  conducted  using  classi¬ 
cal  small  deflection  theory,  and  solutions  were  obtained  only  for  special 
classes  of  plates  which  resulted  in  considerable  simplification  of  the 
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anedysis.  Recent  studies  have  obtained  solutions  to  more  genered  plate 
problems  with  transverse  shear  and  rotary  inertia  effects. 

In  addition  to  the  reports  and  articles,  some  excellent  books  have 
considered  laminated  composite  plates.  Dietz  (2U)  in  19^9  published  a  good 
review  of  much  of  the  early  analytical  work  in  laminated  composites  and 
gave  a  thorough  listing  of  references.  Ashton,  Hedpin,  euad  Petit  (25) 
in  1969  published  an  excellent  primer  for  those  new  in  the  field.  They 
give  a  good  development  and  discussion  of  the  governing  equations.  Ashton 
and  Whitney  (26)  in  1970  published  a  good  review  of  the  progress  in  lami¬ 
nated  plate  theory  and  the  various  solutions  and  solution  techniques. 

1.1.1  Formulation  of  Equations 

Two  early  investigators,  Reissner  and  Stavsky  (l),  analyzed  a  plate 
composed  of  two  orthotropic  lamlnas  of  equal  thickness  aligned  so  that  the 
principal  material  axes  of  the  two  laminas  are  rotated  at  an  angle  of 
+0  and  -9  with  the  plate  axes  (later  referred  to  as  €Uigle-ply  plates). 

They  found  that  even  for  a  small  deflection  theory  analysis,  there  exists 
a  coupling  between  the  transverse  bending  and  inplane  stretching  equations 
that  does  not  exist  for  orthotropic  plates.  They  showed  that  the  coupling 
effect  exists  in  the  boundary  conditions  as  well  as  in  the  governing 
equations. 

Reissner  smd  Stavsky  (l)  formulated  the  governing  equations  as  two 
Uth  order  partial  differential  equations  in  terms  of  the  Airy  stress 
function  and  the  transverse  displacement.  Due  to  the  coupling  of  the 
txo  equations,  they  must  be  solved  simultaneously,  and  Inplane,  as  well  as 
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transverse  boundary  conditions,  must  be  specified  at  each  boundary  to 
obtain  a  particular  plate  solution.  Stavsky  (2)  formulated  the  govern¬ 
ing  equations  as  an  eighth  order  partial  differential  equation  in  terms  of 
the  Airy  stress  function  and  later  (3)  in  the  form  of  three  simultaneous 
partial  differential  equations  in  the  displacements  u,  v,  and  w.  Al¬ 
though  each  of  the  formulations  have  advantages  in  certain  problems, 
they  present  about  the  same  order  of  difficulty  in  obtaining  a  solution. 

Later  contributions  to  the  laminated  plate  formulations  have  been 
made  by  several  authors.  Tsai  and  Azzi  (U)  added  thermal  stresses  to 
the  formiilation  of  laminated  plate  equations.  Whitney  and  Leissa  (5), 
and  Yang,  Norris,  and  Stavsky  (6)  included  roteory  inertia  terms  in  the 
formulation  and  the  latter  authors  also  included  transverse  shear  which 
becomes  important  for  thick  plates.  The  formulation  by  Yeuig,  Norris,  and 
Stavsky  results  in  five  coupled  partial  differential  equations  in  terms  of 
the  displacements  u,  v,  and  w  and  two  rotations  of  the  normals  to  the 
midplane. 

1.1.2  Solutions  for  Symmetric  Plates 

Although  available  formvilations  of  the  governing  equations  and 
boundary  conditions  are  applicable  for  genereO.  laminated  plates,  most 
of  the  solutions  have  been  obtained  for  special  classes  of  plates  which 
result  in  considerable  simplification  to  the  governing  equations. 

Stavsky  (2)  in  his  early  work  showed  that  for  laminated  plates  symmetric 
about  the  geometric  midplane  both  in  properties  and  fiber  orientation, 
the  transverse  bending  and  inplane  stretching  eqiiations  and  the  boundary 
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conditions  become  uncoupled  and  can  be  solved  Independently.  However, 
he  showed  that  even  for  symmetric  laminated  plates,  the  governing  equations 
and  boundary  conditions  still  contain  some  cross-stiffness  terms  due  to 
the  laminas  which  result  in  the  governing  equations  being  different 
from  those  of  an  orthotropic  analysis.  He  further  showed  that  only  for 
very  special  kinds  of  symmetric  plates  do  the  cross-stiffness  terms  vanish 
and  the  governing  equations  become  identical  with  the  orthotropic  equa¬ 
tions.  As  a  result  of  these  simplifications,  symmetric  plates  have  been 
studied  extensively  in  the  literattire. 

The  Rayleigh-Rltz  method  of  solution  with  beam  mode  shape  functions 
has  been  used  effectively  by  several  investigators  to  obtain  deflection, 
vibration,  and  buckling  solutions.  Ashton  emd  Waddoups  (7)  investigated 
symmetric  laminated  plates  and  presented  analytical  results  for  uniaxially 
loaded  plates  with  the  loaded  edges  clamped  and  the  unloaded  edges  free. 
They  also  presented  calculated  frequencies  and  mode  shapes  for  canti¬ 
levered  plates  which  showed  good  agreement  with  experiment.  Ashton  (8) 
extended  the  analysis  to  include  nonuniform  cross  section  auid  material 
properties  and  presented  results  for  tapered  plates  under  unaxial 
loading  with  simple  supported  and  clamped  boundary  conditions.  Ashton 
(9)  further  extended  the  anedysis  to  include  elastically  restrained 
boundary  conditions  and  presented  some  deflection  results  with  re¬ 
strained  boTindary  conditions. 

Later  Srinivas  and  Rao  (10)  obtained  closed  form  solutions  for  the 
free  vibration  and  buckling  of  simply  supported  symmetric  laminated 
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thick  plates  using  linear  small  deflection  theory.  They  ccmipared  their 
results  with  thin  plate  theory  and  showed  that  the  error  associated 
with  thin  plate  theory  increases  with  plate  thickness  and  the  modular 
ratio  between  the  lamina.  A  similar  analysis  by  Whitney  (ll)  showed 
the  same  trends.  Thus  the  definition  of  "thin"  for  laminated  plates 
must  take  into  account  modular  ratio  as  well  as  plate  thickness. 

1.1.3  Solutions  for  Angle-ply  and  Cross-ply  Plates 

Angle-ply  plates  are  unsymmetric  and  consist  of  an  even  number  of 
layers  having  the  same  thickness  and  elastic  properties,  and  the  ortho- 
tropic  ajces  of  symmetry  for  each  ply  are  alternately  orientea  .  angles 
of  +9  and  -9  to  the  plate  axes  (See  fig.  1.).  Cross-ply  plates  are  un¬ 
symmetric  and  consist  of  an  even  number  of  layers  all  of  the  same  thick¬ 
ness  and  elastic  properties  with  the  orthotropic  axes  of  symmetry  in 
each  ply  alternately  oriented  at  0  degrees  and  90  degrees  to  the  plate 
euces  (See  fig.  2j.  Whitney  and  Leissa  (5)  showed  that  for  angle-ply 
plates  and  cross-ply  plates,  the  governing  equations  and  boundary  condi¬ 
tions  are  considerably  simplified  but  are  still  coupled  and  must  be  solved 
simiiltaneously.  As  a  result  of  the  simplifications,  angle-ply  and  cross- 
ply  laminated  plates  have  been  studied  by  severeLL  investigators  to  deter¬ 
mine  the  effects  of  bending-extensioneJ.  coupling. 

Several  solution  procedures  have  been  developed  which  are  effective 
in  solving  the  coupled  governing  equations  and  boundary  conditions. 

Whitney  and  Leissa  (5)  and  Whitney  (12)  used  a  Fovirier  series  technique 
to  obtain  solutions  for  deflection,  free  vibration,  and  buckling  of  both 
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angle-ply  and  cross-ply  plates  for  slnply  supported  boundary  conditions 
which  allow  inplane  displacements  normal  and  tangent  to  the  boundary  for 
cross-ply  and  angle-ply  plates  respectively.  Bert  emd  Mayberry  (13)  used 
a  Raylelgh-Ritz  approach  with  beam  mode  shape  functions  to  obtain  free 
vibration  resvilts  for  cross-ply  and  angle-ply  plates  with  clamped 
boundary  conditions.  Whitney  (lU)  used  a  Galerkin  procedure  to  obtain 
solutions  for  the  shear  buckling  of  cross-ply  plates  with  simply  suppor¬ 
ted  boundary  conditions. 

In  all  cases  cited,  the  bendlng-extensional  coupling  significantly 
affected  the  results.  Bendlng-extensional  coupling  has  the  overall 
effect  of  reducing  the  plate  stiffness  and  thus  increasing  the  static 
deflection  and  reducing  the  natural  frequencies  and  buckling  loads.  The 
effect  of  coupling  is  increased  as  the  number  of  laminas  are  decreased 
and  as  the  degree  of  anisotropy  between  lamina  is  increased. 

Fortier  (15)  used  a  Rayleigh-Ritz  solution  procedure  to  investigate 
the  effects  that  various  types  of  inplane  boundary  restraints  have  on 
the  behavior  of  angle-ply  emd  cross-ply  plates  with  small  initial  curva- 
tvore.  He  also  considered  the  effects  of  transverse  shear  and  inplane 
loads  on  the  natural  frequencies.  He  found  that  inplane  boundary  condi¬ 
tions,  initial  curvature,  transverse  shear,  and  Inplane  loads  all  have 
a  significant  effect  on  the  behavior  of  unsymmetric  plates.  However, 
Whitney  (l6)  considered  transverse  shear  for  cylindrical  bending  of 
symmetric  and  cross-ply  plates  and  showed  that  treuisverse  shear  has  less 
effect  on  the  deflection  of  cross-ply  plates  than  symmetric  plates. 
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1.1.4  Approximate  Solutions 

Due  to  the  difficulty  of  obtaining  general  plate  solutions,  approxi¬ 
mate  solutions  have  emerged  also.  A  "reduced  bending  stiffness"  method 
was  proposed  by  Ashton  (20)  in  vhlch  the  bending-extenslonal  coupling 
is  accounted  for  in  an  approximate  wfky.  Approximate  solutions  were  ob¬ 
tained  by  reducing  the  bending  stiffness  of  the  plate  by  an  amount  deter¬ 
mined  by  the  coupling  terms  and  then  neglecting  the  coupling  effect. 
Ashton  used  the  method  to  compare  maximxss  deflections  of  a  cross-ply  euid 
an  angle-ply  plate  with  measured  results  cmd  obtained  good  agreement. 
Whitney  (21)  also  used  the  method  to  compare  results  vith  those  calcu¬ 
lated  using  the  Fourier  series  technique  and  shoved  good  agreement  for 
both  angle-ply  and  cross-ply  plates  for  simply  supported  boundary  condi¬ 
tions,  but  for  clamped  boundary  conditions,  relatively  poor  agreement  was 
obtained  for  certain  euigle-ply  cases. 

A  more  refined  anedysls  was  presented  by  Srlnivas  (22)  which  con¬ 
siders  transverse  sheeo*  deformations  emd  rotary  inertia  effects.  By 
assuming  piecewise  linear  variations  of  the  displacements  u  and  v 
and  constant  values  of  w  across  the  plate  thickness,  the  problem  be¬ 
comes  two  dimensional.  Trigonometric  series  solutions  were  obtained  for 
the  approximate  two  dimensional  problem  for  plates  with  singly  supported 
boundazy  conditions,  and  the  results  showed  good  agreement  with  exact 
theory. 


1.1.5  Dynamic  Stability 

A  limited  nimiber  of  investigators  have  considered  the  dynamic 
stability  of  laminated  plates.  Smirnov  (17)  considered  the  flutter  of 
an  infinite  sandwich  plate  subjected  to  cylindrical  bending  in  a  gas 
stream  and  obtained  an  exact  solution.  He  later  extended  the  analysis 
(l8)  to  include  semi-infinite  sandwich  plates  with  various  boundary  con¬ 
ditions.  He  used  linear  piston  theory  aerodynamics  and  obtained  solu¬ 
tions  for  both  clamped  and  simply  supported  semi-infinite  plates. 

Librescu  and  Badoiu  (19)  and  Ramkumar  (23)  analyzed  the  flutter  of 
flat,  symmetrically  laminated,  siii^>ly  supported  plates  with  inplcme 
normal  loads.  They  used  linear  piston  theory  dynamics  with  aerodynsunic 
damping  and  the  Raylei^-Ritz  solution  procedure  with  beam  node  shape 
functions.  Flutter  boundaries  are  presented  as  a  function  of  the 
orientation  of  the  principal  axis  of  the  lamina  for  various  inplane 
loads  and  aerodynamic  dancing  coefficients.  However,  it  should  be 
noted  that  all  the  flutter  boundaries  presented  in  reference  19  were 
obtedned  using  only  two  terms  in  both  the  streamwise  and  cross-stream 
direction  and  the  results  presented  in  reference  23  were  obtedned  using 
ten  terms  in  the  streamwise  direction  but  only  two  terms  in  the  cross¬ 
stream  direction.  Thus,  in  both  cases,  the  results  presented  may  not 
have  been  completely  converged. 

1.2  Statement  of  Problem 

Althou^  much  work  has  been  done  in  the  analysis  of  composite 
laminated  plates,  as  discussed  in  the  literature  survey,  panel  flutter 
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stands  out  as  having  received  very  little  attention.  Panel  flutter  has 
long  been  recognized  as  a  problem  for  the  design  of  conventional  panels 
sxibjected  to  supersonic  flow.  This  is  evident  by  the  large  number  of 
reported  panel  flutter  investigations  some  of  which  are  listed  in  refer¬ 
ences  27-U2.  The  panel  flutter  work  that  has  been  done  for  laminated 
plates  is  so  limited  that  it  is  of  little  value  to  the  designer.  It  is 
evident  that  to  design  efficient  laminated  conqposite  panels  for  super¬ 
sonic  application,  a  better  understanding  of  their  flutter  characteris¬ 
tics  is  needed. 

Another  area  that  hats  received  little  attention  is  the  solution  of 
general  laminated  plate  problems  wher^  the  plates  may  have  any  number 
of  layers  stacked  in  an  arbitrary  sequence  with  the  fibers  in  each  layer 
rotated  at  an  arbitrary  angle  to  the  plate  axis.  For  general  laminated 
plates,  the  conq)lete  bending  and  extensional  governing  equations  and 
boundary  conditions  are  coupled  and  must  be  solved  simultaneously  to 
obtain  plate  solutions.  Considerable  difficulty  is  encountered  in 
solving  the  coupled  equations  and  boundary  conditions.  Thus,  althou^ 
the  governing  equations  and  boundary  conditions  have  been  formulated  for 
general  laminated  plates,  no  general  solutions  were  found  in  the  litera- 
tvu:e  siiTvey.  All  the  solutions  presented  were  for  specially  laminated 
(symmetric,  angle-ply,  or  cross-ply)  plates  whose  special  geometry 
results  in  sin^lifi cations  to  the  governing  equations  eind  boundary  con¬ 
ditions  which  make  it  considerably  easier  to  obtain  solutions.  However, 
since  in  the  design  of  laminated  plates  it  is  not  always  practical  or 
desirable  to  use  specially  laminated  plates,  an  analysis  is  needed  which 


obtains  solutions  for  general  laminated  plates. 

The  purpose  of  this  investigation  is  to  obtain  the  flutter  charac¬ 
teristics  of  flat,  general  laminated  plates  using  approximate  methods. 

The  analysis  will  be  based  upon  small  deflection  theory  but  will  consider 
the  general  coupled  governing  equations  and  boundary  conditions  for 
simply  supported  plates.  Thus,  an  approximate  solution  procedure  will 
be  developed  for  con^jletely  general  plates  which  may  have  any  nimber  of 
layers  stacked  in  an  arbitrary  sequence  and  with  the  fibers  in  each 
lamina  rotated  at  an  arbitrary  angle  to  the  plate  axes.  The  plate  may 
be  subjected  to  uniform  inplane  normal  and  shear  loads. 

An  extended  Galerkln  method  will  be  used  to  obtain  approximate 
solutions  to  the  governing  equations  and  boundary  conditions.  The 
extended  Galerkin  method  wiJl  be  used  since  it  provides  a  straight¬ 
forward  solution  procedure  for  nonconservative  problems  (See  ref.  Ua) 
using  simple  series  to  describe  the  displacements.  The  aerodynamic 
loading  on  the  panel  will  be  e^sumed  to  be  given  by  linear  piston  theory 
£Uid  the  flow  may  be  at  an  arbitrary  cross-flow  angle.  Piston  theory 
aerodynamics  will  be  used  because  it  gives  simple  expressions  for  the 
aerodynamic  loading  and  has  been  shown  to  be  applicable  for  Mach  numbers 
greater  than  1.6  (See  ref.  28  and  29i}. 

Since  symmetric  and  angle-ply  plates  have  been  the  subject  of  many 
investigations,  approximate  solutions  will  be  obtedned  for  typical 
laminated  plates  from  each  class,  emd  a  limited  parametric  study  will  be 
conducted  to  determine  the  effects  on  the  flutter  boundaries.  The  param¬ 
eters  to  be  studied  will  include  the  number  and  orientation  of  the  plies. 
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length-width  ratio,  inplane  shear  and  normal  loads,  and  cross-flow 
angles.  The  flutter  boundaries  from  this  analysis  will  be  coII^>ared  with 
those  calculated  using  clsissical  orthotropic  plate  theory  and  the  reduced 
bending  stiffness  method. 

Since  general  laminated  plates  have  not  been  investigated  prior  to 
this  time,  approximate  solutions  of  their  governing  eqimtions  will  be  of 
special  interest.  Thus,  flutter  boxmdaries  will  be  calculated  for  some 
typical  general  laminated  plates  and  the  results  will  be  compared  with 
those  obtained  for  specially  constructed  symmetric  and  angle-ply  plates. 
Flutter  boundaries  will  also  be  calculated  for  some  composite  stiffened 
aluminum  plates  that  do  not  have  a  symmetric  or  angle-ply  construction. 

Althoiagh  the  approximate  procedure  will  be  developed  for  the  purpose 
of  obtaining  flutter  solutions,  the  analysis  will  edso  have  the  capability 
of  giving  natural  vibration  frequencies,  and  inplane  normal  and  shear 
static  buckling  loads.  Since  no  numerical  results  are  available  in  the 
literature  for  general  laminated  plates,  a  significant  contribution  to 
the  literature  could  be  made  by  using  the  natural  vibration  and  static 
buckling  capabilities  of  the  analysis.  However,  only  a  limited  number 
of  natural  frequency  calcvilations  will  be  made  to  compare  with  pviblished 
results  to  verify  the  solution  procedure. 


Chapter  2 


DEVELOPMENT  OF  GOVERNING  EOUATIONS 

In  order  to  develop  the  governing  equations  for  this  investigation, 
several  assumptions  are  made  as  follows: 

1.  The  plate  is  constructed  of  flat,  uniform  thickness  layers  of 
oirthotropic  sheets  bonded  together.  The  direction  of  principad  stiffness 
of  the  indivldvial  layers  do  not  in  general  coincide  with  the  plate  edges. 

2.  The  plate  is  thin;  l.e.,  the  thickness  is  much  smaller  than  the 
other  physical  dimensions. 

3.  The  displacements  are  small  compared  to  the  thickness. 

4.  Each  lamina  obeys  Hooke's  lav. 

5.  The  Kirchhoff  hypothesis  is  used;  i.e.,  normals  to  the  midplane 
of  the  undeformed  plate  remain  straight  and  normal  to  the  midplane  during 
deformation. 

6.  Transverse  shear  and  nomal  strains  ai*e  negligible. 

7.  Body  and  rotary  inertia  forces  cure  negligible. 

8.  The  plate  is  of  constant  thickness. 

These  assumptions  give  rise  to  the  conclusions  that  y  ,  y  ,  T  ,  T  ,  and 

xz  yz  xz 

a  are  negligible  which  is  the  case  for  eui  approximate  state  of  plane 
z 

stress. 

The  coordinate  system  used  to  identify  the  plate  and  lamina 
geometry  is  shown  in  figure  3.  The  distances  to  the  Indlvidued  laminas 
eire  meastired  from  the  geometrical  mldpleuie  of  the  plate.  The  positive 
directions  for  stresses  are  shown  in  figure  4. 


2.1  Lamina  Stress-Strain  Equations 


From  elementary  strength  of  materials,  Hooke's  law  for  an  orthotro¬ 
pic  lamina  in  a  state  of  plane  stress  is  given  as  follows  (See  ref. 


(2.1) 


where  the  stresses  ®i»®2*^12*  strains  "^2.2 

referred  to  the  direction  of  principal  stiffness  and  the  subscript  I 
refers  to  a  particular  lamina.  This  relationship,  written  for  another 
(x-y-z)  system  of  axes  (See  figure  5.)  where  the  x-y  axes  are  rotated  at 

an  angle  0  with  respect  to  the  1-2  axes,  is  given  as  follows: 

5.1  S.2 

5.2  ^22 

where  the  Q^j  matrix  is  a  transformed  matrix  with  the  elastic  proper¬ 
ties  along  the  prlnclpel  axes  of  the  lamina  rotated  to  the  x-y-z 
(plate)  axes  system.  The  matrix  has  non-zero  values  for  all  terms 
when  the  princlpcLl  axes  of  the  lamina  do  not  coincide  with  the  plate 
axes  systan  (for  example,  a  fibrous  composite  with  the  fibers  rotated  at 
an  angle  with  the  plate  axes). 
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From  elementary  strength  of  materieLLs  considerations,  the  strain  at 
any  point  in  a  lamina  undergoing  deformation  can  be  expressed  in  terms 
of  the  deformation  of  the  geometrical  midplane  of  the  lamina.  For 
small  deflections  euad  a  state  of  plane  stress,  the  strain  at  any  point 
in  a  lamina  z  distance  from  the  midplane  is  given  in  terms  of  the 
midplane  strain  and  curvature  as  follows: 
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where  e  °,  e  ®,  and  y  °  are  the  midplane  strains  and  <  ,K  ,  and  K 

X  ’  y  *  'xy  x’  y  xy 

are  the  midplane  ciurvatures.  The  midplane  strains  and  cvirvatures  de¬ 
fined  in  terms  of  the  displacements  are  needed  in  the  later  development 
of  the  equations  and  are  given  as  follows: 
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Using  equations  (2.2)  and  (2.3) •  the  stress  state  at  any  point  in  a 
lamina  may  be  written  In  terms  of  the  midplane  strain  and  curvature  as: 
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»  [Q]^  +  z  [Q]j^  [<]j^ 


where  z  is  the  dlsteunce  from  the  midplane  to  the  point. 
2.2  Lamina  and  Plate  Constitutive  Equations 


(2.6) 


U  I 


(2.7) 


Since  for  a  plate  cooqosed  of  several  laminas,  the  stresses, 
strains,  and  displacements  will  be  different  for  each  lamina,  it  is 
convenient  to  define  some  equivalent  system  of  forces  and  moments  which 
will  be  considered  to  be  applied  to  the  midplane  of  the  plate .  Within 
the  plane  stress  assumption,  such  a  system  will  be  defined  for  a  lamina 
in  terms  of  three  stress  and  three  moment  resultants.  The  stress  and  moment 

resultants  for  a  lamina  axe  defined  as: 
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and  are  statically  equlvedent  to  the  actual  stresses  on  the  lamina  and 
may  be  considered  to  be  applied  at  the  midplane  of  the  lamina.  The 
positive  directions  for  the  stress  and  moment  resultants  on  the  plate 
are  shown  in  figure  6.  Using  equations  (2.6),  (2.8)  and  (2.9)  the 
stress  and  moment  resultants  can  be  related  to  the  strain  and  curvature 
of  the  midplane.  The  resulting  relations  are  known  as  the  lamina 
constitutive  equations  and  are  ^xven  as  follows: 
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The  plate  constitutive  equations  will  be  developed  from  equations  (2.10) 
and  I (2.11)  by  summing  up  the  indivlduea  terms  for  each  lamina. 

Thus, 
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Since  e**  emd  ic  are  functions  of  x  and  y  only  and  Q  Is  a  function  of 
lamina  properties,  they  can  be  taken  outside  the  integred..  Thus, 
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or  after  integrating  these  may  be  written  as: 
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where 
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Writing  the  constitutive  equations  for  laminated  plates  in  matrix 
notation,  they  become: 

[N]  -  [A]  [e®]  +  [B]  [k]  (2.19) 

[M]  =  [B]  [e®]  +  [D]  [k]  (2.20) 

These  equations  indicate  that  for  a  general  lamixiated  plate  either  stress 
or  moment  resultants  produce  both  strains  and  curvatures.  Thus  there  is 
a  coupling  between  the  stress  and  moment  resultants  through  the  matrix  [B]. 
Since  [B]  is  an  even  function  of  h  (see  eq.  (2.18),),  it  is  zero  and  hence 
the  equations  are  \incoupled  for  classical  orthotropic  plates  cmd  for 
laminated  plates  that  are  symmetric  about  the  plate  midplane. 

Equations  (2.19)  and  (2.20)  will  be  rewritten  in  a  more  useful  form 
for  the  following  development  by  solving  equation  (2.19)  for  the  midplane 

strains,  and  rewriting  equation  (2.20)  in  terms  of  the  stress  resultants. 
Equations  (2.19)  and  (2.20)  become: 


(a) 


(b) 

(2.18) 


(c) 
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[e®]  -  [A*]  [N]  [B*]  [k1  (2.21) 

[M]  =  In]  +  [D*]  [k]  (2.22) 

where 

[A*]  »  [A"^] 

[B*]  =  -[A"^]  [B] 

[D«]  =  [D]  -  [B]  [A"^]  [B]  (2.23) 

2.3  Plate  Governing  Equations 

Due  to  the  coupling  between  moaent  and  stress  resxiltants  noted  In 
equations  (2.19)  and  (2.20),  two  differential  equations  are  required  to 
describe  the  behavior  of  the  plate  under  various  loaid  conditions.  The 
first  equation  will  be  developed  using  plate  equations  of  motion  and  the 
second  equation  will  be  developed  from  the  compatibility  conditions  that 
must  be  satisfied  between  the  strains. 

2.3.1  Equations  of  Motion 

The  equations  of  motion  for  a  thin  laminated  plate  cure  identical  to 
those  of  homogenous  plate  theory.  For  a  plate  subjected  to  Inplane  normal 
and  shear  loads  and  loads  perpendicular  to  the  plate  surface,  the  equations 
of  motion  are  given  as  follows  (See  ref.  26^  where  the  inertia  terms  In  the 
X  and  y  directions  are  considered  to  be  negligible. 


0 


(2.2U) 
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an  dN 


ay  ax 


(2.25) 


aSi  aSi 


a^ 

*  ^ 


-  p  -  N, 


"  ax2 


!I  ^  -  2  N 

y  3y2  xy 


a2  a2 

+  V  ^ 


(2.26) 


The  inplane  loads  N  .  N  ,  and  N  may  be  separated  into  the  uni- 

X  y  xy 

form  externally  applied  loads  and  the  loads  induced  by  the  plate  deflections 
as  follows: 


N  *  »  ♦ 

X  XX 


N  =  N  +  N' 

y  y  y 


N  =  N  +  N’ 
xy  xy  xy 


(2.27) 


Substituting  the  expressions  for  N  ,  R  ,  and  N  given  by  equation  (2.27) 

X  y  xy 

into  equations  (2.24),  (2.25) t  >nd  (2.26)  results  in  the  following 
equations  of  motion: 


aN’ 

aN' 

—SL 

3y 

aN* 

aN'  ■ 

3y 

ax 

(2.28) 


(2.29) 
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3x‘ 


3y' 


3^ 


2  *  ^  3x3y 


P-f. 


a2 

3  V 


(2.30) 


.2  .2 

_  N  li  +  V  IJi 
V  2  '  2 

y  3y^  3t 


*  3x^ 


The  higher  order  terms 
have  been  neglected  In  equation  (2.30). 

Defining  a  stress  function  F  such  that: 


.2  ^2  «2 
/»!  3  V  »•  3  V  j  wi  3  V  » 

(»..  O  .  O  .  “d  g^gy)  , 


3y 


N' 

X 


2 


— 5  (a) 


3y 


N*  = 

y 


^  (b) 

3x^ 


(2.31) 


»• 

xy 


3^F 

'5^ 


(c) 


then  equations  (2.28)  and  (2.29)  are  satisfied  Identically  by  the  function. 
The  function  F  Is  referred  to  In  the  literature  (See  ref.  26. )  as  the 
Airy  stress  function.  Using  the  definitions  of  F  given  by  equations 
(2.3l)t  definitions  of  curvature  given  by  equation  (2.3) »  and  the 
expression  for  the  moments  given  by  equation  (2.22)  In  equation  (2.30), 
the  following  expression  Is  obtained: 


(2.32) 
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a 
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32v 
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3x2  ■ 

3*  2  D* 

^22  3y2  “  ^  “26  ?x5y 
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_  a2  ^2  ^2 

T^  w  g  w  K  3  w  o  w  a  V  . 
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»  9x2 


or 

*  3V  .  •  • 

®21  ®26  -  ®6l 


)  a^F 


3x^3y 
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3x^3y 


♦  2  (D*  ♦  2  D*  )  aI»-2  ♦  U  D*  ^*d' 
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^2  ~2  kc 

p*H  H*a»  l-L 

*  9x2  y  X, 


3x3y- 
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(2.33) 


If  p  is  taken  as  the  aerodynaniic  loading  due  to  air  flov  over  the 
plate  surface,  then  equation  (2.33)  is  one  of  the  equations  needed  to 
solve  for  the  flutter  of  flat  laminated  plates.  Aerodynamic  loads  pre¬ 
dicted  by  piston  theory  have  often  been  used  in  the  literature  (See  ref. 
28  -  U2O  to  obtedn  flutter  solutions.  Piston  theory  gives  a  relatively 
simple  expression  for  p  and  has  been  shown  to  be  applicable  for  Mach 

niasbers  greater  than  1.6  (See  ref.  28  and  29.).  For  supersonic  flow 
over  one  side  of  a  flat  plate  at  an  angle  of  A  with  the  x-axls  as 
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sbovn  In  figure  7,  piston  theory  gives  the  foUovlng  expression  for  the 
aerodynamic  loads. 

P  ‘  cos  A  +  ^  sin  Aj  {2.3k) 

Substituting  the  expression  for  p  given  by  equation  {2.3k)  into 
equation  (2.33)»  the  following  governing  equation  is  obtained  for  flutter 
of  flat  laminated  plates. 


B 


•  3^ 

/  *  •  \ 

♦ 

3**F  /  * 

3  *  <®U  * 

3x^ay 

4  -  ^  4> 

2  (2.35) 

ax  ay 
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■  '2  «16 
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”  ®62^  3x  3y^ 

3x^3y 

,  •  »  »  3**w 

"  2  (D,2  -  2  Dgg) 


^  3‘*w 

^2  *'‘^26  — 


3x‘3y 


3  +  ^22  ^ 
3x3y^  3y^ 
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cos  A  + 


—  3^w  3^ 

W  9x3y  ^  3^ 


This  expression  is  a  Uth  order  partial  differential  equation  which  governs 
the  flutter  behavior  of  flat  laminated  plates  but  is  in  terms  of  two  un¬ 
knowns  w  euid  F.  Thus  an  additioncLl  equation  relating  w  and  F  is 
needed  before  solutions  to  flutter  problems  may  be  obtained. 


ro|< 


2.3.2  Compat ibillty  Equation 


The  additional  equation  needed  to  solve  for  the  plate  flutter 
behavior  is  obtained  from  the  coa^atibility  condition  that  must  be 
satisfied  between  the  midplane  strains.  The  appropriate  equation  is 
obtained  by  eliminating  the  midplane  displacements  from  the  strain- 
displacement  relations  given  in  equation  (2.U).  Differentiating 
(eq.  2.U(a))  twice  with  respect  to  y  and  e°  (eq.  2.U(b))  twice  with 
respect  to  x  and  adding,  the  following  expression  is  obtained: 


aV  ^  aV 

axay^  ax^ay 


(a) 


(2.36) 


but  from  eq.  (2.U(c)): 


aV  ^  aV 

axay^  ax^ay 

therefore ; 


(b) 


(2.37) 


Equation  (2.37)  is  referred  to  as  the  compatibility  equation  and  will  be 
used  to  develop  the  second  governing  equation  in  terms  of  w  and  F. 


Using  the  expressions  for  strains  given  'by  equation  (2. 21)  and  the 
force  resultants  given  by  equations  (2.27)  and  (2.31)  and  the 
curvatures  ,•  given  by  equation  (2.5)  in  equation  (2.37),  give  rise  to  the 
following  exjtresslon: 
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?r-  *  ®22  -  2  Bgg)  3,V 
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(2.39) 


This  is  the  second  governing  equation  needed  to  solve  for  the  plate  be¬ 
havior.  It  is  also  a  fourth  order  partial  differentieQ.  equation  in  tenas 
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of  F  and  w  and  for  flutter  problems  must  be  solved  simultaneously  with 
equation  (2.35).  In  order  to  solve  the  equations  for  a  specific  problem 
the  boundary  conditions  must  be  specified. 

2.k  Boundary  Conditions 

Proper  bovtndary  conditions  which  gucurantee  unique  solutions  to  the 
two  governing  equations  must  be  sx>eclfied.  It  has  been  found  (See 
reference  26.)  that  the  necessary  boundary  conditions  are  those  of 
classical  homogeneous  plate  theory  plus  those  of  an  inplane  elasticity 
problem.  Thus  due  to  the  coupling,  the  usual  deflection,  slope,  moment,  or 
shear  boundary  conditions  used  in  classifical  plate  theory  do  not  give 
unique  solutions  but  inplane  boundary  conditions  must  also  be  sx>ecified. 

The  resxxlting  boundary  conditions  require  one  member  of  each  pair  of  the 
following  quantities  to  be  specified  along  the  boundaries: 


w 

or 

3M  ^ 
nt 

3t 

♦  N  |2  . 

n  3n 

N  ^ 
ns  3s 

(a) 

(2.U0) 

3w 

3n 

or 

M 

n 

(b) 

u® 

n 

or 

N 

n 

(c) 

or 

N  4. 

nt 

(d) 

where  n 

and 

t  8u:e 

used 

to  denote 

coordinates 

normcd. 

and  tangential  to  the 

plate  boundary  respectively. 
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For  this  analysis,  the  boundary  conditions  chosen  represent  those 
for  a  simply  supported  plate  with  no  inplane  edge  restraints  eind  are 
given  as  follovs: 


V  (0,y)  »  w  (a,y)  «  w(x,0)  -  w(x,b)  »»  o  (a) 

\  \  ^a»y)  *  0  (^) 

M  (x.O)  »  M  (K,b)  ■  0  (c) 

O'  Jr 

(O.y)  -  (a.y)  *  0  (d) 

N*  (x.O)  =  N*  (x.b)  »  0  (e) 

Or  jr 


(O.y)  »  (a.y)  -  N^(x.O)  »  N^(x.b)  »  0  (f) 


(2.U1) 


These  equations  will  be  rewritten  in  terms  of  w  and  F  using  the  expression 
for  moments  given  by  equation  (2.22)  and  the  definitions  of  strain,  cur¬ 
vature,  and  F  given  by  eq^ationa  (2.h),  (2.5),  and  (2.31)  respectively. 
Thus  the  boundary  conditions  become: 

w  (O.y)  =  w  (a.y)  =  w  (x.O)  ■  w(x,b)  »  0  (2.42) 
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(2.44) 
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2 

d  F 

— P  «  0  at  X  ■  0,a 

3y 


0  at  y  =  0,b 


(2.U5) 


(2.U6) 


3^ 

3x3y 


0  at  X  «  0,a 
y  ■  0,b 


(2.1*7) 


The  governing  equations  (2.35)  and  (2.39)  along  with  the  boundary  con¬ 
ditions  (2.1*2)  through  (2.1*7)  will  be  solved  by  using  approximate  techniques 
in  the  next  chapter. 


Chapter  3 

APPROXIMATE  SOLUTIONS  OF  THE  FLUTTER  EQUATIONS 

The  two  governing  fourth  order  partial  differential  equations  (2.35) 
and  (2.39)  edong  with  boundary  conditions  given  by  equations  (2.U2) 
through  (2. 1+7)  will  be  used  to  obtain  approximate  flutter  solutions.  The 
analysis  will  be  developed  for  the  flutter  of  general  laminated  plates  where 
arbitrary  stacking  sequence  and  orientation  of  the  lamina  fibers  is 
permissible.  Since  certain  special  classes  of  plates  have  been  studied 
almo?5t  exclusively  in  the  literature,  the  simplifications  to  the 
general  solution  due  to  the  special  constructions  will  be  discussed.  The 
special  classes  of  plates  considered  are  symmetric  laminated  plates  and 
angle-ply  plates.  Symmetric  plates  may  be  composed  of  any  even  number  of 
lamina  of  arbitrary  thickness  and  orientation  of  the  fiber  directions  as 
long  as  for  every  lamina  above  the  plate  midplane  there  is  an  identical 
(in  thickness,  material  properties,  cuid  orientation)  lamina  equal 
distance  below  the  plate  mldplane.  Angle-ply  plates  are  less  general  than 
symmetric  plates  and  consist  of  an  even  number  of  laminas  all  of  the  same 
thickness  and  elastic  properties  with  the  orthotropic  axes  of  symmetry 
in  each  lamina  alternately  oriented  at  +Q  and  -6  to  the  plate  axis. 

Although  the  analysis  and  solution  procedure  will  be  developed  basically 
to  study  panel  flutter,  the  analysis  is  general  and  the  resulting  computer 
program  can  also  be  used  to  determine  static  buckling  loads  Md  natural 
vibration  frequencies. 


31 


32 


3.1  General  Laminated  Plates 

An  extended  Galerkln  method  will  he  used  to  obtain  approximate 
flutter  solutions  since  it  provides  a  strai^t-foivard  solxition  proce- 
dvire  for  nonconservative  problems  (See  ref.  U2.)  using  a  simple  series 
to  describe  the  assumed  displacements.  The  extended  Galerkin  method 
admits  solutions  only  when  the  ass\amed  deflections  satisfy  the  geometric 
(deflection  and  slope)  boundary  conditions.  However,  the  number  of 
terms  required  in  the  solution  to  obtsdn  converged  results  are  probably 
reduced  if  the  assumed  deflections  satisfy  some  or  all  of  the  natural 
(force  or  natural  constraints)  boundsiry  conditions.  Since  functions  are 
not  available  which  satisfy  all  the  boundary  conditions  given  by  equa¬ 
tions  (2.k2)  through  (2.k7),  the  capability  to  account  for  the  natural 
boundary  conditions  makes  the  extended  Gederkin  method  particularly 
suited  to  this  analysis. 

The  extended  Galerkin  method  of  solution  is  illustrated  by  the 
following  equation  written  in  terms  of  the  virtual  work  for  displace¬ 
ments  w. 


b  a  a  b  a 

J  f  (governing  equation)6v  dxdy  *  f f^(w,x,y)  6w  jdx  +  f 
0*1)  0  0  0 


t  a.  "b  a 

*  J  6w.j  dy  +  j  fj^(w,x,y)  6  j  dy  =  0  (3.1) 

0  0  0  0 


dx 
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The  terms  under  the  double  and  single  Integrals  are  usually  individueLIly 
set  equal  to  zero  emd  referred  to  as  the  governing  equation  and  boundary 
conditions,  respectively.  For  the  boundary  conditions  to  be  satisfied, 
either  the  functions  f(w,x,y)  or  the  variation  of  the  deflections  or  the 
slopes  must  be  zero  at  the  appropriate  bovindaries.  In  usual  applications 
of  the  Galerkin  method,  a  series  of  functions  is  chosen  which  satisfy  all 
the  boundary  conditions  and  the  coefficients  of  the  functions  are  deter¬ 
mined  so  that  the  double  integral  term  is  zero  or  equivalently  the  coeffi¬ 
cients  are  determined  so  that  each  term  of  the  series  is  orthognal  to  the 
exact  solution  (See  ref.  U3O.  When  all  the  natural  boundary  conditions 
are  not  satisfied,  the  extended  Galerkin  method  requires  that  the  coeffi¬ 
cients  be  determined  using  the  unsatisfied  single  integral  terms  as  well 
as  the  double  integrad  terms  in  equation  3.1« 

For  the  flutter  of  general  laminated  plates,  a  series  of  functions 
for  both  w  and  F  are  needed  which  exactly  satisfy  the  geometrical 
boundary  conditions  given  by  equation  (2.1*2)  and  which  satisfy  as  closely 
as  possible  the  natural  boundary  conditions  given  by  equations  (2.1+3) 
through  (2.1*7).  The  functions  assumed  for  this  analysis  are  given  as 
follows : 


m=l  n=l 


(3.3) 


These  functions  exactly  satisfy  the  geometrical  boundary  conditions  and 
the  natiural  boundary  conditions  given  by  equation  (2.U7)  but  the  re¬ 
maining  naturaJ.  boundary  condltioBs  are  not  completely  satisfied.  The 
unsatisfied  natural  boundary  condition  terms  are  given  as  follows: 
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The  boundary  condition  terms  that  must  be  Included  in  the  extended 
Galerkin  method  written  in  terms  of  the  virtual  work  are  given  as  follows 
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Ny: 


■) 


dx 


(3.8) 


where  the  displacements  u°  and  are  given  as  follows  and  were  obtained  by 
substituting  the  definition  of  curvature  given  by  equation  (2.5)  and  the 
expression  for  strain  given  by  equation  (2.22)  into  equation  (2.^)  and 
integrating. 
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Using  the  assumed  functions  given  by  equations  (3.2)  and  (3.3)  and  applying 
the  extended  GeuLerlLln  method,  the  orthogonality  relations  for  the  governing 
equations  (3.35)  and  (3.39)  with  the  boundaiy  condition  terms  given  by 
equations  (3.5)  through  (3.10)  added  may  be  written  as  follows: 
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o  x**l  s*X  I  I 
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Integrating  and  rearranging  the  terms,  the  following  equations  are 
obtained; 


39 


vhere 


(3.16) 


(3.17) 


Equations  (3.lU)  and  (3.15)  results  in  2()(xN)  equations  in  terns  of  the 

two  unknowns  and  In  order  to  get  the  equations  in  a  form 

that  can  be  solved*  H _  must  be  solved  for  in  terms  of  C  .  This  can 

ra  rs 

be  done  using  equation  (3.15)  and  writing  it  in  the  following  matrix 
form; 


(3.18) 


where  r  and  s  are  summed  from  1  to  M  and  N  respectively  and 

coefficient  matrices  for  H  and  C  res- 
mnrs  mnrs  mn  mn 

pectively.  To  obtain  H  as  a  function  of  C  ,  premultiply  each  side 

ITS  rs 

of  equation  (3.18)  by  the  Inverse  of  fb  1  as  follows: 
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Equation  (3.20)  results  in  a  series  solution  for  H  in  terms  of  the 

mn 

unknown  coefficients  C  . 

rs 
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Using  equation  (3.20)  to  eliminate  from  equation  (3.lU),  MxN 

homogenous  equations  are  obtained  in  terms  of  the  unknown  C  .  The 

irS 

resulting  equations  may  be  written  in  the  following  form: 


where  the  c  are  the  coefficients  of  C _ .  Equation  (3.21)  is  recog- 

mn  mn 

nized  as  the  well  known  characteristic  equation  and  has  a  non-trivial 
solution  obtedned  by  setting  the  (Jeterminant  of  the  coefficient  matrix 
to  zero.  This  results  in  an  algebraic  equation  of  order  MxN  in  terms 
of  w.  The  roots  of  the  algebraic  equation  are  the  characteristic  veLlues 
or  eigenvalves  and  represent  the  plate  nondimensional  vibration  fre¬ 
quencies  . 

Flutter  solutions  have  been  obtedned  using  a  digital  computer  to 
solve  equation  (3.21)  for  its  elgenvedves .  Since  the  coefficients. 


c . . ,  are  a  function  of  the  flutter  parameter,  X,  the  plate  vibration 

i  J 

frequencies  are  also  a  function  of  X.  The  criteria  used  to  define  the 


point  of  flutter  is  the  lowest  value  of  X  which  results  in  two  of  the 


frequencies  (w^^j)  coalescing  and  thus  having  a  negative  imaginary 
value.  This  criteria. was  selected  because  a  negative  imaginary  fre- 
qviency  results  in  the  assumed  deflection  function  (See  eqviation  (3.2).) 


Ul 


becoming  unbounded  since  it  is  multiplied  by  e  .  Similar  criteria 
has  been  used  extensively  in  the  literature  (See  ref.  39-)  to  define  the 
point  of  flutter. 

3.2  Symmetric  Laminated  Plates 


Plates  that  are  laminated  such  that  they  are  symmetric  about  the 

midplane  represent  the  largest  class  of  laminated  plates  of  practical 

use.  Their  symmetric  construction  results  in  warp-free  structural 

elements  which  are  desirable  in  most  applications.  Also,  the  governing 

differential  equations  emd  the  boundary  conditions  are  considerably 

simpler  than  for  general  laminated  plates.  Thus,  symmetric  plates  have 

been  investigated  more  than  euiy  other  class  of  laminated  plates. 

As  was  pointed  out  in  section  (2.2)  and  as  shown  by  equation 

(2.18(b)),  B  is  an  even  function  of  the  leunina  thickness  and  for 
mn 

symmetric  plates  all  the  terms  vanish.  Thus,  the  flutter  equations 

(equations  (2.35)  and  (2.39))  are  considerably  shortened  and  more  impor¬ 
tantly  become  uncoupled.  Then  only  equation  (2.35)  needs  to  be  solved 
to  obtain  a  flutter  solution.  For  a  symmetrically  laminated  plate,  the 


governing  flutter  equation  becomes; 


3x  3y 
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terns  since  D 


where  the  asterisks  have  been  left  off  the  D 
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for  symmetric  laminated  plates. 


Although  equation  (3.22)  is  considerably  sitopler  than  the  flutter 
equations  for  general  laminated  plates,  it  differs  from  the  governing 
equation  for  flutter  of  an  orthotropic  plate  by  the  additional  and 

terms.  These  terms  appear  in  the  governing  equation  as  a  result  of 
the  fibers  not  being  aligned  with  the  plate  axes. 

Since  for  symmetric  plates  the  governing  equations  become  uncoupled, 
the  boundazy  conditions  that  must  be  satisfied  for  a  flutter  solution 
are  given  by  equations  (2.U2)  throu^  (2.U1*)  and  are  rewritten  els 


follows  with  the  B  terms  removed: 

mn 
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(2.23) 

(2.2U) 


(2.25) 


These  equations  differ  from  the  boundary  conditions  for  a  simply 
supported  orthotropic  plate  by  the  ^16  and  terms . 

To  obtain  a  standard  Oalerkln  solution  for  symnetric  plates.  It  Is 
necessary  to  assiuie  a  function  for  v  that  satisfies  all  the  boundary 
conditions  (3.23)  throuj^  (3.25).  However,  the  and  Dgg  terms 

make  it  difficult  to  .find  satisfactory  functions.  Therefore,  the  same 
deflection  function  used  in  section  (3*1)  for  genercLL  laminated  plates 
and  given  by  equation  (3.2)  will  be  assumed  and  the  extended  Oalerkln 


1*3 


method  lased  to  account  for  the  unsatisfied  natural  boundary  conditions. 
The  boundary  condition  terms  that  must  be  accounted  for  by  the  extended 
Galerkin  method  are  given  as  follows: 
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Using  the  assumed  functions  given  by  equation  (3.2)  and  applying  the 
extended  Galerkin  method  to  the  governing  equation  (3.22)  with  the 
and  Dgg  boundary  condition  terms  included,  results  in  the  following 


eqmtion; 
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Equation  (3.29)  represents  a  system  of  (MxN)  linear  homogenous  equations 

in  terms  of  the  unknown  C  and  is  recognized  as  the  eigenvalue  equa- 

mn 

tion.  Flutter  solutions  are  obtained  for  equation  (3.29)  in  the  same 
manner  as  discussed  in  section  (3.1)  for  general  laminated  plates. 

3 • 3  Angle-ply  Plates 


Angle-ply  plates  by  definition  are  constructed  so  as  to  have  an  even 
number  of  layers  all  of  the  same  thickness  and  elastic  properties  an^ 
with  the  orthotropic  axis  of  symmetry  in  each  ply  alternately  oriented 
at  +6  and  -0  to  the  plate  axis  (See  Fig.  2.).  Although  angle-ply 
plates  do  not  represent  as  large  or  as  in5x>rtant  a  group  as  the  symmetric 
plates,  they  are  inqpoirtant  from  an  analytic  viewpoint  for  obtaining  a 
better  understanding  of  general  laminated  plates.  Because  of  the  special 
geometry  of  angle-ply  plates,  the  governing  equations  are  considerably 
simplified  but  are  still  coupled.  Thus,  angle-ply  plates  retain  many  of 
the  characteristics  of  general  laminated  plates  but  relatively  simple 
solutions  ceui  be  obtedned.  These  characteristics  have  resulted  in  the 
angle-ply  plates  being  studied  extensively  in  the  literature  to  obtain  a 
better  understanding  of  general  laminated  plates. 

For  angle-ply  plates,  it  can  be  shown  (See  ref.  5.)  that  this 
special  construction  results  in  some  of  the  coefficients  being  zero  as 
follows : 
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Setting  these  terms  to  zero  in  the  flutter  equations  (2.35)  and  (2.39), 
results  in  the  following  governing  equations. 
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These  equations  axe  coupled  and  must  he  solved  simultaneously. 

Since  the  governing  equations  for  angle-ply  plates  are  coupled,  all 
the  boundary  conditions  given  by  equations  (2.U2)  through  (2.U7)  must  be 
satisfied.  However  with  some  of  the  coefficients  being  zero  as  indicated 
by  equations  (3.29),  the  boundary  conditions  are  simplified  as  follows: 


w(0,y)  =  w(a,y)  =  w(x,0)  =  w{x,b)  =  0 
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If  the  functions  for  w  and  F  used  for  the  general  laminated 
plates  (section  3.1)  and  given  by  equations  (3.2)  emd  (3.3)  are  used, 
then  boundary  conditions  (3.33),  (3.3^),  (3.35),  and  (3.38)  are  satis¬ 
fied.  Boundary  conditions  (3.36)  and  (3.37)  are  not  satisfied  directly 
but  vill  be  shovm  to  be  satisfied  indirectly  by  the  following  develop¬ 
ment.  Since  equations  (3.36)  and  (3.37)  are  not  satisfied,  the  terms 
that  must  be  included  in  the  extended  Galerkin  method  are  given  by 
equations  (3.7)  and  (3.8)  and  are  repeated  as  follows: 


(3.7) 


(3.8) 


Using  equations  (3.9) »  (3.10),  (3.29)  and  (3.30),  the  displacements  u° 
6uid  v°  are  given  as  follows: 
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Using  the  assumed  functions  given  by  equations  (3.2)  suid  (3.3)  in  equa¬ 
tions  (3.39)  and  (3.^0)  results  in  zero  values  for  u°  and  v°  on  the 
boundaries  x  =  0  and  a,  and  y  =  0  and  b,  respectively.  Thus,  all 
the  boundary  conditions  are  satisfied  without  using  the  extended 
Galerkin  method. 

Substituting  the  expressions  for  F  and  w  given  by  equations 
(3.2)  and  (3.3)  into  equations  (3.31)  and  (3.32)  results  in  the  follow¬ 
ing  equations : 
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Equation  (3A5)  results  in  MxM  homogenous  equations  in  terms  of 
and  is  recognized  as  the  eigenvalue  equation.  Flutter  solutions  are 
obtained  for  equation  (3.^5)  in  the  same  manner  as  discussed  in  section 
3.1  for  general  and  symmetric  laminated  plates. 

3.U  Approximate  Solutions 

An  approximate  theory  (See  ref,  20.)  referred  to  as  the  reduced 
bending  stiffness  theory  will  be  coii5)ared  with  the  present  analysis.  In 
the  approximate  theory,  the  plate  bending  stiffness  is  reduced  by  an 
amount  depending  on  the  co\Q>llng  between  the  governing  equations  eind 
then  the  coupling  is  neglected  when  solving  the  equations.  The  reduced 
bending  stiffness  is  the  same  as  that  defined  by  equation  (2.18(c)). 
Thus,  the  governing  equation  for  the  reduced  bending  stiffness  theory 
is  obtained  from  equation  (2.35)  by  neglecting  the  coupling 
is  given  as  follows:  • 
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Since  the  coupling  is  neglected,  only  boundary  conditions  given  by  equa¬ 
tions  (2.42)  through  (2.44)  with  [B  ]  “  0  need  to  be  satisfied.  The 

nin 

boundary  conditions  are  given  as  follows: 

w(0,y)  =  w(a,y)  =  w(x,0)  =  w(x,b)  «  0  (3.47) 
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The  resxolting  governing  eqiiation  (3.46)  and  boundary  conditions 
(3.47)  through  (3.49)  are  identical  to  those  used  for  symmetric  lamin¬ 
ated  plates  (See  equations  (2.22)  through  (2.25).)  except  that  [D*]  is 
used  instead  of  [O].  Thus,  the  approximate  theory  reduces  the  solution 
for  general  laminated  plates  to  that  of  a  symmetric  plate  with  the  bend¬ 
ing  stiffness  reduced.  For  angle-ply  plates;  the  and  D*26 

(See  equation  (3.29).)  ^.^e  zero  and  the  approximate  governing  equation 
(3.46)  and  boundary  conditions  (3.47)  throu^  (3.49)  are  identical  to 
those  for  classical  plate  theory.  Thus,  the  reduced  bending  stiffness 
theory  can  be  used  to  obtain  solutions  for  angle-ply  plates  from 
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piibllshed  classical  plate  theory  results.  Since  for  symmetric  plates 
[B^]  =  0  and  [D*]  =  [d],  the  reduced  bending  stiffness  theory  and  the 
present  analysis  are  identical. 

The  advantages  of  using  the  approximate  theory  are  evident  for 
angle-ply  plates  since  results  can  be  obtained  using  veil  knovn  published 
solutions.  For  symmetric  and  genereil  laminated  plates,  the  advantages 
are  not  so  clear  since  a  symmetric  analysis  similar  to  the  one  discussed 
in  section  3.2  is  necessary  even  to  obtain  approximate  solutions.  How¬ 
ever,  prior  to  the  present  analysis  the  reduced  bending  stiffness  theory 
was  the  only  solution  procedure  available  for  general  laminated  plates. 


Chapter  U 

RESULTS  AND  DISCUSSIONS 


The  equations  goverhlng  the  behavior  of  general  laminated  plates 
have  been  programed  for  approximate  solution  on  a  digital  cooqputer. 
Althou^  the  computer  programs  were  developed  largely  to -obtain  flutter 
results,  natural  -vibration  frequencies  and  static  buckling  loads  may 
readily  be  obtained  using  the  programs.  Special  analysis  and  programs 
were  de-veloped  to  obtain  solutions  for  symmetric  and  angle-ply  plates. 
However,  the  ad-vantages  of  using  the  special  progreuns  over  the  one 
developed  for  the  general  laminated  plates  were  minimal  and  in  all  cases 
gave  identical  results. 

Flutter  boundaries  have  been  calculated  for  symmetric,  angle-ply, 
and  general  laminated  plates.  The  plate  nomenclature  smd  geometry  are 
shown  in  figure  8  for  symmetric  plates,  in  figure  9  for  euagle-ply  plates, 
and  in  figure  10  for  general  laminated  plates.  Flutter  boundaries  were 
also  obtEdned  for  an  aluminum  plate  with  one,  two  and  four  layers  of 
coiiq>osite  material  applied  to  one  or  both  sides  of  the  plate  as  shewn  in 
figvu'e  11.  The  material  properties  of  the  individual  lamina  used  in 
making  the  cedculations  are  typical  of  those  for  boron-epoacy  and  glass- 
epoxy  matericLLs.  These  material  properties  are  given  in  table  1. 

Material  properties  for  the  composite  stiffened  aliuninum  plate  are  gi-ven 
in  table  II. 
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U.l  Convergence  of  Results 

When  using  a  series  solution.  It  is  In^ortant  to  determine  that 
sufficient  terms  are  used  in  the  analysis  to  obtain  converged  results. 
However,  using  more  terms  than  necesseury  results  in  considerably  longer 
computation  times.  The  usual  procedure  for  determining  convergence  is 
to  start  with  a  small  number  of  terms  find  increase  the  number  until  the 
solution  does  not  change  as  more  terms  sire  included.  For  solutions 
obtained  herein,  a  converged  result  is  one  which  changes  no  more  than 
1  percent  when  more  terms  axe  included  in  the  flutter  determinant.  This 
method  of  determining  convergence  is  demonstrated  in  figure  12  for  a 
four-ply,  symmetric,  square  panel  where  the  flutter  parameter  is  shown 
as  a  function  of  the  number  of  terms  in  the  x-direction.  Ctirves  are 
shown  for  0=0°  and  60°,  and  the  symbols  represent  the  different 
numbers  of  terms  used  in  the  y-direction.  For  0  =  0°,  fo\ir  terms  in 
the  x-direction  and  only  one  tern  in  the  y-direction  are  sufficient  for 
converged  results.  For  0  =  60°,  approximately  six  terms  in  both  the 
X-  and  y-directions  are  needed  to  obtain  results  converged  within  1  per¬ 
cent.  Note  that  for  0  =  60°,  the  converged  solution  is  approached  from 
hi£^er  values  of  X  for  m  >  k  and  n  >  2.  This  suggests  that  the 
cross-stiffness  terms  have  an  adverse  effect  on  the  flutter  solution,  a 
phenomena  that  will  be  discussed  in  a  later  section. 

The  modes  that  coeLLesce  to  produce  flutter  for  0  =  0°  and  0  =  60° 
are  shown  in  figures  •13(a)  and  13(b),  respectively,  where  X  is  shown 
as  a  function  of  the  non-dimensional  frequency  squared.  For  both  cases. 


the  (1,1)  and  (2,1)  modes  coalesce  and  resxilt  in  flutter.  The  different 
nunhers  of  terms  required  to  obtain  converged  results  for  the  two  cases 
are  due  to  the  ^16  and  Dgg  cross-stiffnesstenas.  For  0  =0°,  the  cross¬ 
stiffness  terms  are  zero,  and  the  plate  is  an  orthotropic  plate  for 
which  it  has  been  shown  (See  ref.  UO.)  that  only  the  first  mode  in  the 
cross-stream  direction  has  an  effect  on  the  flutter  solutions.  For 
9  =  60°,  the  cross-stiffness  terms  couple  the  modes  in  the  x-  and 
y-directions,  and  consequently,  all  modes  effect  the  flutter  solutions. 
An  attempt  has  been  made  throughout  the  remainder  of  this  study  to  in¬ 
sure  that  the  values  presented  are  converged  within  approximately  1 
percent . 


k.2  Comparison  with  Literature 


In  order  to  verify,  as  well  as  possible,  the  accuracy  of  the  present 
analysis,  the  euialysis  will  be  \ised  to  obtain  flutter  data  comparable 
with  that  presented  in  reference  19  for  square  symmetric  plates.  The 
material  properties  for  the  plate  considered  in  reference  19  are  given 
in  table  III.  A  comparison  of  the  flutter  boundaries  is  shown  in  figure 
lU  \rtiere  X  is  shown  as  a  function  of  the  orientation  angle  of  the  out¬ 
side  lamina.  The  dashed  curves  were  obtained  from  reference  19,  and  the 
solid  curves  were  obtained  using  the  present  analysis.  Both  sets  of 
curves  were  calculated  using  only  two  terms  in  the  x-  and  y-directlons . 
The  agreement  between  the  two  analyses  is  very  good  except  for  0'  ®  0° 
and  0'  =  90°  and  for  0  near  90°  where  large  discrepancies  do  exist. 
Since  the  present  analysis  uses  a  Galerkin  method  of  solution  and  the 
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einalysis  of  reference  19  uses  a  Rayleig^-Ritz  method,  the  discrepancy  is 
attributed  to  the  difference  between  the  two  methods  when  an  insufficient 
number  of  terms  are  used  in  the  analysis  to  obtain  converged  results. 

For  converged  results,  both  einalyses  should  give  the  same  values.  Con¬ 
verged  flutter  boundaries  obtained  using  the  present  analysis  are  shown 
in  figure  15.  Although  it  is  evident  that  the  two-term  solutions  are 
not  converged,  there  is  still  good  agreement  between  the  trends  shown  by 
the  two-term  and  the  converged  analysis.  For  6=0°  or  0  =  90°  auid 
for  6*=  0°  or  0'  =  90°,  the  present  analysis  should  give  the  same  re¬ 
sults  as  an  orthotropic  analysis  since  the  cross-stiffness  terms  are 
zero.  This  is  shown  in  figure  15  where  the  symbols  represent  orthotropic 
results  obtained  from  reference  31. 

Although  no  flutter  solutions  were  found  in  the  literature  for 
8uigle-ply  or  general  laminated  plates,  some  comparisons  of  natural  fre¬ 
quencies  for  angle-ply  and  cross-ply  plates  enre  possible.  The  four 
lowest  natural  frequencies  calculated  using  the  present  analysis  are 
compared  in  tables  IV  eind  V  with  those  calculated  in  reference  15  vising 
a  Rayleigh-Ritz  method  for  square  angle-ply  and  cross-ply  plates,  res¬ 
pectively.  The  natural  frequencies  calculated  by  the  present  analysis 
for  simple  supports  should  compare  directly  with  those  referred  to  in 
reference  15  as  having  hinge-free-tangential  and  hinge-free-normal 
boundary  conditions,  respectively,  for  the  angle-ply  and  cross-ply 
plates.  For  angle-ply  plates,  natural  frequencies  are  given  for  the 
principal  material  directions  rotated  at  0°,  10°,  20°,  30°,  and  U0° 

with  respect  to  the  plate  axes.  The  agreement  between  the  two  analyses 
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is  excellent. 

For  cross-ply  plates,  natural  frequencies  calculated  using  the  two 
analyses  are  compared  (See  table  V.),  for  values  of  E^/Egg  =  1,  10,  20, 
30,  and  Uo.  The  agreement  between  the  two  analyses  is  excellent  except 
for  the  higher  frequencies  and  high  values  of  the  ratio  ®2.l'^^22' 
ever,  the  differences  are  too  small  to  be  of  practical  significance. 

Since  a  special  analysis  was  not  made  for  cross-ply  plates,  the  analysis 
for  general  laminated  plates  was  \ised  to  obtain  the  natureil  frequencies. 
Thus,  the  assumed  mode  shape  functions  do  not  satisfy  the  boundary  con¬ 
ditions  exactly,  but  they  are  satisfied  throu^  the  extended  Galerkin 
procedure  as  discussed  in  section  (3.1).  The  good  eigreement  between  the 
natural  frequencies  indicates  the  validity  of  using  the  Galerkin  method 
to  account  for  the  boundary  conditions  not  satisfied  by  the  assumed  mode 
functions. 

U.3  Symmetric  Plates 

Flutter  boxindaries  for  symmetric  plates  with  lamina  properties 
typical  of  boron-epoxy  material  or  a  glass-ejjoxy  material  are  presented 
in  this  section  where  A  is  shown  sis  a  function  of  the  angle  that  the 
fibers  make  with  the  x-axis.  The  effects  of  inplane  normal  euid  shear 
loads  and  cross-flow  on  the  flutter  of  a  square  boron-epoxy  plate  com¬ 
posed  of  k  laminas  are  also  shown.  The  values  of  A  in  each  case  are 
based  on  which  is  the  vedue  of  with  all  the  fibers  aligned 


with  the  x-axis. 
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4. 3.1  Flutter  Boundaries 

Flutter  boundaries  for  square,  symmetric,  glass-epoxy  and  boron- 
epoxy  plates  are  presented  in  figure  l6(a),  and  l6(b),  respectively,  for 
plates  composed  of  two,  four,  and  six  laminae.  For  square  symmetric 
plates,  the  highest  flutter  stability  is  obtained  with  the  fibers  aligned 
with  the  x-axis;  rotating  the  fibers  away  from  the  x-axis  results  in  a 
continuoxis  reduction  in  flutter  stability  for  values  of  9  up  to  90°. 
Flutter  boundaries  calculated  vising  cleissical  orthotropic  plate  theory 
which  neglects  the  cross-stiffhess  terms  ^16  and  are  seen  to  be 

inaccurate  and  nonconservative.  Thus,  inclusion  of  the  cross-stiffness 
terms  in  the  euialysis  has  a  destabilizing  effect  on  flutter.  Increeislng 
the  number  of  laminas  composing  the  plate  results  in  inqxroved  flutter 
stability  and  flutter  boundaries  that  are  closer  to  those  given  by 
classicaLL  plate  theory.  A  comparison  of  the  flutter  boundaries  for  the 
glass-epoxy  (See  fig.  16(a).)  and  boron-epoxy  (See  fig.  l6(b).)  plates 
shows  that  increasing  the  orthotropy  (for  glass-epoxy  material 


=  3,  euid  for  boron-epoxy  materieLL  °^  laminas  causes 

the  discrepancy  between  classical  plate  theory  euid  the  present  analysis 
to  be  more  pronovinced.  Also,  for  the  more  highly  orthotropic  material 
properties,  orienting  the  fibers  transverse  to  the  flow  has  a  larger 
destabilizing  effect  on  flutter. 

Flutter  bovindarles  similar  to  those  presented  in  figure  l6  for  a 
square  plate  are  presented  in  figure  17  for  a  boron-epoxy  plate  with  a 
length-width  ratio  of  2.0.  These  flutter  boundaries  show  the  solutions 
based  upon  classical  plate  theory  to  be  hi^ly  nonconservative  due  to 
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the  neglect  of  the  cross-stiffness  terms.  For  the  a/b  =  2.0  plate, 
rotation  of  the  fibers  away  from  the  x-axis  results  in  an  increase  in 
the  flutter  stability  to  maximum  values  of  X  for  6  of  approximately 
30°  and  40°,  respectively,  for  plates  composed  of  two  and  four  laminas. 
Thvis,  in  designing  composite  laminated  plates,  proper  fiber  orientation 
is  important  to  improved  flutter  stability  for  the  plate. 

U.3.2  Effect  of  Inplane  Normal  and  Shear  Loads 

The  effect  of  the  flutter  boundaries  of  inplane  normal  and  shear 
loads  for  a  square,  symmetric,  boron-epoxy  panel  composed  of  four 
laminas  is  shown  in  figures  l8  and  19  where  X  is  presented  as  a  func¬ 
tion  of  the  inplane  load.  Flutter  boundaries  are  shown  for  values  of 
0=0°,  15°,  and  30°.  The  circle  symbols  indicate  the  point  of  buckling 
with  air  flow. 

Inpleine  normal  loads  (See  fig.  18.)  resiilt  in  approximately  a  linear 
reduction  in  flutter  stability  for  loads  iq)  to  the  point  of  buckling. 

The  effects  of  the  cross-stiffness  terms  and  rotation  of  the  fibers  away 
from  the  x-axis  both  result  in  a  reduction  in  the  flutter  stability  and 
buckling  loads.  Thus,  classical  plate  theory,  which  neglects  the  cross¬ 
stiffness  terms,  is  nonconservative  in  its  predictions  of  both  flutter 
stability  and  buckling  loads. 

Small  inplane  shear  loads  (See  fig.  19.)  have  a  stabilizing  effect 
on  flutter  for  0  =  15°  and  0  =  30°,  but  leurger  loads  decrease  the 
flutter  stability.  The  increase  in  flutter  stability  with  small  inplane 
shear  loads  is  contrary  to  results  obtained  for  isotropic  plates  (See 
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ref.  Ul.)  and  is  due  to  a  change  in  the  modes  that  coalesce  to  produce 

flutter.  This  is  shown  in  fig\ire  20  where  coalescence  of  the  frequencies 

for  6  =  15°  are  shown  for  values  of  inplane  shesu"  loads  =  0 . ,  .  5 , 

and  1.0.  For  K  =0  (See  fig.  20(a).)  and  K  =  1.0  (See  fig.  20(c).), 
xy  xy 

the  two  lowest  frequencies  which  correspond  to  the  (1,1)  and  (l,2)  inodes 

coalesce  and  result  in  flutter.  However,  for  a  small  inplane  shear  load 

K  of  0.5  (See  fig.  20(b).),  the  lowest  frequency  and  one  of  the  higher 

xy 

frequencies,  which  correspond  to  the  (l,l)  and  (2,l)  modes,  coalesce  to 

produce  flutter.  For  K  greater  thEui  1.0,  the  (l,l)  and  (1,2)  modes 

xy 

continue  to  coalesce  to  produce  flutter  but  at  lower  values  of  X  (See 
fig.  19.).  The  cross-stiffness  terms  and  rotating  the  fibers  away  from 
the  x-axis  both  result  in  increased  buckling  loads  and  in  improved 
flutter  stability  for  most  values  of  inplane  shear  loads.  This  is  con¬ 
trary  to  results  noted  for  inplane  normal  loads.  It  should  be  noted 
that  at  the  point  of  buckling,  the  flutter  values  of  A  for  inplane 
shear  loads  (See  fig.  19.)  are  considerably  lower  than  those  of  inplane 
normal  loads  (See  fig.  18. ) .  This  is  in  contrast  with  results  presented 
in  reference  itl  for  isotropic  plates. 

U.3.3  Effect  of  Cross-flow 

The  effects  of  cross- flow  on  the  flutter  boundaries  for  a  square, 
symmetric,  boron-epoxy  plate  coinpcsed  of  four  laminas  are  shown  in 
figure  21.  Flutter  boundaries  are  presented  for  0=0°,  15°  and  30° 
where  A  is  shown  as  a  function  of  the  cross-flow  angle.  For  9  =  15° 
and  6  =  30°,  flow  at  small  cross-flow  angles  may  have  a  substantial 
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stabilizing  effect  on  flutter.  For  6  =  15°,  a  30  percent  increeise  in 
X  may  be  obtained  for  a  10°  cross-flow  angle  whereas  larger  cross-flow 
angles  resvilt  in  a  reduction  in  flutter  stability.  Similar  results  are 
obtained  for  0  =  30°  where  a  45  percent  increase  in  X  i^  obtained 
for  a  cross- flow  angle  of  20°.  This  increase  in  the  flutter  stability 
with  cross-flow  angle  can  be  explained  as  follows:  Cross-flow  is  de¬ 
stabilizing  for  an  orthotropic  plate  (See  ref.  40.).  Orienting  the 
fibers  at  an  angle  with  the  plate  axis  results  in  the  maximum  bending 
stiffness  and  th\is  the  maximim  flutter  resistance  occurring  in  that 
direction.  These  factors  tend  to  covmteract  each  other  and  result  in 
the  maximxim  flutter  stability  occurring  at  an  angle  between  the  x-axis 
and  the  fiber  directions.  Note  that  the  maximum  value  of  X  for  the 
cases  Shown  are  obtained  at  an  angle  approximately  two-thirds  of  the 
angle  A  the  fibers  make  with  the  x-axis. 

For  plates  with  the  fibers  aligned  with  the  x-axLs  (6  =  0°),  cross- 
flow  is  destabilizing,  but  at  sma.ll  cross-flow  singles  they  have  higher 
values  of  X  than  plates  with  the  fibers  rotated  at  0  =  15°  and  0  = 
30°.  Classical  plate  theory  does  not  show  any  beneficisLi  effects  of 
cross-flow  and  thus,  underpredicts  the  flutter  stability  for  large  cross- 
flow  angles.  For  laminated  cos^osite  plates,  cross-flow  may  be  either 
stabilizing  or  destabilizing  depending  on  the  flow  angle  and  the  orienta¬ 
tion  of  the  fibers.  Thus,  in  designing  composite  plates  that  will  experi¬ 
ence  cross-flow,  an  optimum  orientation  of  the  fibers  may  considerably 
iii5)rove  the  flutter  stability  of  the  plate. 
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U.U  Angle-ply  Plates 

Flutter  boundaries  calculated  vising  the  present  analysis,  classical 
plate  theory,  and  reduced  bending  stiffness  theory  are  presented  in  this 
section  for  typical  angle-ply  plates.  The  effects  of  inplane  normal  and 
shear  loads  and  cross-flow  on  the  flutter  of  a  square  boron-epoxy  plate 
composed  of  four  laminas  are  also  shown.  There  is  a  bending-extensional 
coupling  between  the  governing  equations  for  an  angle-ply  plate  which  is 
neglected  by  classical  plate  theory.  Since  the  cross-stiffness  terms 
discussed  in  section  U.3  for  symmetric  plates  are  zero  for  angle-ply 
plates,  the  only  difference  between  classical  plate  theory  and  the  pre¬ 
sent  analysis  is  due  to  bending-extensional  coupling.  The  reduced  bend¬ 
ing  stiffness  theory  accounts  for  the  coupling  in  an  approximate  way  by 
reducing  the  bending  stiffness  of  the  plate  by  an  amount  determined  from 
the  coupling  terms  and  then  neglecting  the  coupling  in  solving  the  equa¬ 
tions.  For  0=0°  or  0  =  90°,  the  coupling  terms  are  zero,  and  the 
present  analysis  and  the  reduced  bending  stiffness  theory  become  identi¬ 
cal  with  classical  plate'  theory. 

k.k.l  Flutter  Boundaries 

Flutter  boundaries  for  square,  angle-ply  plates  with  properties 
typical  of  a  glass-epoxy  and  a  boron-epo:^  material  eure  presented  in 
figure  22(a)  and  22(b),  respectively,  where  X  is  shown  as  a  function 
of  0.  Bending-extensional  coupling  has  a  large  destabilizing  effect  on 
the  plate  composed  of  only  two  laminas,  but  as  the  nvimber  of  laminas 
increase,  the  coupl-.xig  effect  becomes  smaller,  and  the  boundaries 
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approach  those  given  by  classical  plate  theory.  The  reduced  bending 
stiffness  theory  gives  flutter  boundaries  that  are  in  good  agreement 
with  the  actual  boundaries  even  for  the  plate  composed  of  only  two 
laminas.  The  small  errors  associated  with  using  the  reduced  bending 
stiffness  theory,  however,  are  not  always  conservative.  Comparing  the 
flutter  boundaries  for  the  glass-epoxy  (See  fig.  22(a).)  and  the  boron- 
epoxy  (See  fig.  22(b).)  plates  shows  that  Increasing  the  orthotropy  of 
the  laminas  Increases  the  destabilizing  effect  of  bending-extensional 
coupling.  The  good  agreement  between  the  present  analysis  and  the  re¬ 
duced  bending  stiffness  analysis  does  not  deteriorate  as  the  lamina 
orthotropy  is  increased. 

Flutter  boundaries  similar  to  those  presented  in  figure  22  for  a 
square  plate  are  presented  in  figure  23  for  a  boron-epoxy  plate  with  a 
length-width  ratio  of  2.0.  The  flutter  boundaries  show  larger  destabil¬ 
izing  effects  due  to  bending-extensional  coupling  than  fovind  for  the 
sqviare  plates,  but  increasing  the  number  of  laminas  conposing  the  plate 
reduces  the  destabilizing  effectj.  For  plates  with  four  or  more  laminas, 
orienting  the  fibers  at  an  angle  with  the  x-axis  may  result  in  large 
inprovements  in  the  flutter  stability  with  the  highest  flutter  values 
being  obtained  at  9  -  U5°.  For  plates  with  two  laminas,  rotating  the 
fibers  away  from  the  x-axis  results  in  little  or  no  improvement  in  the 
flutter  stability.  In  each  case,  orienting  the  fibers  transverse  to  the 
x-axis  results  in  the  most  unstable  condition.  Good  agreement  is  obtained 
between  the  present  analysis  and  the  reduced  bending  stiffness  theory 
for  all  plates,  but  the  largest  discrepancy  is  obtained  for  the  plate 
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composed  of  two  laminas. 

4.4.2  Effect  of  Inplane  Normal  and  Shear  Loads 

The  effects  of  inplane  normal  and  shear  loads  on  the  flutter  of 
square  angle-ply  plates  are  shovn  in  figure  24(a)  euid  24(b),  respectively. 
Flutter  boundaries  are  shown  for  values  of  0  =  0°,  15°,  and  30°.  The 
circle  symbols  indicate  the  point  of  buckling  with  air  flow.  Inplane 
normal  and  shear  loads  result  in  a  shcurp  drop  in  the  flutter  stability 
for  loads  up  to  the  point  of  buckling.  Bending-extensional  coupling  has 
a  destabilizing  effect  on  flutter  and  reduces  the  psoiel  buckling  loads. 
For  plates  with  inplane  normal  loads  (See  fig.  24(a).),  rotating  the 
fibers  away  from  the  x-axis  generally  has  a  destabilizing  effect  on  the 
flutter  and  reduces  the  buckling  loads,  but  for  plates  with  large  in¬ 
plane  shear  loads  (See  fig.  24(b).),  the  opposite  effect  occurs.  How¬ 
ever,  for  small  values  of  inplane  shear,  rotating  the  fibers  away  from 
the  x-axis  may  also  be  destabilizing.  The  reduced  bending  stiffness 
theory  gives  flutter  boundaries  that  are  in  good  agreement  with  the  pre¬ 
sent  analysis  for  each  case  shown.  It  should  be  noted  that  at  the  point 
of  buckling,  the  flutter  values  of  A  for  Inplane  shear  loads  are  con¬ 
siderably  lower  than  those  for  inplane  normal  loads.  This  is  in  agree¬ 
ment  with  results  shown  for  symmetric  laminated  plates. 

4.4.3  Effect  of  Cross-Flow 

The  effect  of  cross-flow  on  the  flutter  boundcurles  for  a  square, 
angle-ply,  boron-epoxy  plate  composed  of  four  laminas  is  shown  in  figure 
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25.  Flutter  boundaries  are  presented  for  6=0°,  15°,  and  30°  where  X 
is  shown  as  a  function  of  the  cross-flow  angle.  In  each  case,  cross- 
flow  and  bending-extensional  coupling  have  a  large  destabilizing  effect 
on  the  flutter  boundary.  However,  rotating  the  fibers  away  from  the  x- 
axis  results  in  cross-flow  having  a  less  destabilizing  effect  on  the 
flutter  boundary.  Thus,  for  angle-ply  plates,  aligning  the  plate  axis 
with  the  flow  results  in  the  most  stable  condition,  but  if  cross-flow  is 
unavoidable,  proper  orientation  of  the  fiber  directions  may  reduce  its 
destabilizing  effects.  For  each  condition  shown,  the  reduced  bending 
stiffness  theory  shews  good  agreement  with  the  present  analysis. 

U.3  General  Laminated  Plates 

The  term  "general  laminated  plates"  as  used  in  this  study  enconH 
passes  plates  composed  of  any  number  of  laminas,  with  arbitrary  thick¬ 
ness,  stacking  sequence,  and  material  properties,  that  satisfy  the 
lineetr  small  deflection  theory  assumptions.  This  definition  covers  an 
infinite  number  of  plates  including  the  symmetric  and  angle-ply  plates 
discussed  in  sections  U.3  and  4.4. 

Flutter  boundaries  will  be  shown  for  two  square,  boron-epoxy  panels 
that  are  conq)osed  of  four  laminais  stacked  in  the  sequence  shown  in  figure 
10  and  designated  as  P-1  and  P-2.  In  cuidition,  flutter  boundeurles  will 
be  calculated  for  six  square,  composite-stiffened  alumintun  plates  con¬ 
figured  as  shown  in  figure  11.  The  flutter  boundaries  calculated  for 
plates  P-1  and  P-2  using  the  present  analysis  will  be  conqiared  with 
those  calculated  using  classical  plate  theory,  and  rec’ 
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stiffness  theory  and  with  those  presented  for  the  8ynB»tric  and  angle- 
ply  plates  in  sections  1;.3  and  h.h.  The  difference  between  the  present 
analysis  and  classical  plate  theory  for  general  laminated  plates  is  due 
to  both  cross-stiffness  terns  and  bending-extensional  covipling. 

U.5.1  Flutter  Boundaries 

Flutter  boundaries  calcvilated  using  the  present  analysis  axe  com¬ 
pared  with  those  calculated  using  classical  plate  theory  and  reduced 
bending  stiffness  theory  in  figure  26(a)  and  26(b),  respectively,  for 
plates  P-1  and  P-2.  The  bending-extensional  coupling  euid  cross-stiffness 
terms  have  a  large  Influence  on  the  flutter  boundary  for  both  plates  and 
render  the  classiceLL  plate  theory  hi^ly  nonconservative.  For  these 
cases,  neglecting  the  coupling  and  cross-stiffness  terms  may  result  in 
flutter  values  that  are  close  to  100  percent  too  hi^^.  However,  flutter 
boundnries  calc\tLated  using  the  reduced  bending  stiffness  theory  show 
good  agreement  with  the  present  analysis.  The  agreement  is  slightly 
better  for  plate  P-1  than  for  '  This  is  expected  since  plate  P-1  is 
more  nearly  symmetric  and  th  .e  coupling  effect  is  smaller  than  for 

jilate  P-2. 

U.S.2  Effect  of  Plate  Construction 

Flutter  boundaries  for  plates  P-1  suid  P-2  are  compared  in  figure  27 
with  those  obtained  for  the  symmetric  and  angle-ply  plates.  The  ciirves 
shown  represent  flutter  boundeurles  for  six  plates,  all  of  which  are  com¬ 
posed  of  the  same  material  and  have  the  same  thickness  but  which  have 
different  flutter  characteristics  as  a  result  of  the  nianber  of  layers 
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and  the  sequence  in  vhich  they  are  laialnated.  Although  the  flutter 
boundaries  for  all  the  plates  follow  the  same  basic  trends  with  6,  the 
angle-ply  plate  with  K  «  U  shows  considerably  more  resistance  to 
flutter  than  any  of  the  other  plate  constructions.  The  symmetric  plate 
with  K  =  U  shows  the  next  best  resistance  to  flutter  followed  by  the 
genereil  laminated  plates  P-1  and  P-2  which  have  similar  flutter  boun¬ 
daries.  The  symmetric  and  angle-ply  plates  with  K  -  2  show  the  lowest 
values  of  X,  which  are  ets  much  as  Uo  percent  lower  than  the  maximum 
values  obtained  for  the  angle-ply  plate  with  K  =  U.  Thus,  improvements 
in  the  flutter  stability  are  obtained  by  increasing  the  number  of  laminas 
and  stacking  them  in  an  angle-ply  sequence. 

The  effects  of  cross-flow  on  the  flutter  boundaries  for  the  general 
laminated  plates  may  be  seen  in  figure  28  where  X  is  plotted  as  a  func¬ 
tion  of  the  cross-flow  angle.  The  flutter  boundaries  for  the  symmetric 
and  angle-ply  plates  with  K  *  U  are  also  shown  for  comparison.  For 
the  general  laminated  plates,  cross-flow  at  small  emgles  slightly  improve 
the  flutter  stability.  For  plate  P-1,  cross-flow  at  angles  greater  them 
20°  and  for  plate  P-2,  cross-flow  at  angles  greater  than  5°  each  resvilts 
in  a  steady  decrease  in  flutter  stability.  Similar  results  are  obtained 
for  the  symmetric  plate  where  large  improvements  in  the  flutter  stability 
are  obtained  for  cross-flow  angles  up  to  10°.  However,  for  the  angle- ply 
plate,  even  a  small  amount  of  cross-flow  results  in  a  reduction  in  flutter 
stability.  Although  it  was  pointed  out  previously  that  for  no  cross-flow, 
the  angle-ply  construction  is  the  most  resistant  to  flutter,  the  presence 
of  cross-flow  may  result  in  the  symmetric  construction  giving  better 
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flutter  stability.  Thus,  in  designing  laminated  plates,  the  stacking 
sequence  may  have  a  significant  effect  on  the  flutter  stability. 

U.5.3  Conqjosite  Stiffened  Aluminum  Plates 

Flutter  boundaries  are  presented  in  figure  29  for  six  square,  com¬ 
posite  stiffened  aluminvua  plates  constructed  as  shown  in  figure  11  and 
with  material  properties  as  given  in  table  II.  Two  of  the  plates  are  of 
symmetric  construction  with  one  or  two  layers  of  conposlte  material 
applied  to  each  side  of  the  CLLuminum.  The  other  four  plates  are  of  vin- 
symmetric  construction  with  one  emd  two  layers  of  coiqposite  material 
applied  to  one  or  both  sides  of  the  aliiminum.  In  each  case,  50  percent 
of  the  plate  mass  is  conqjosed  of  composite  materials,  and  the  only 
difference  between  the  plates  is  the  way  in  which  the  composite  material 
is  applied.  The  flutter  boundaries  for  the  composite  stiffened  plate 
are  coiiq)ared  with  the  flutter  boundary  for  an  equal  mass  aluminum  plate. 
The  values  of  A  shown  in  figure  29  are  based  on  of  the  aluminum 

plate  so  that  the  relative  advantages  of  the  varioxis  constructions  may 
be  seen. 

The  composite  stiffened  aluminm  plates  have  considerably  higher 
flutter  stability  than  the  aluminum  plate.  For  small  values  of  6, 
aluminum  plates  with  lamincis  on  both  sides  of  the  plate  have  values  of 
A  that  are  50  percent  hi^er  than  those  obtained  for  plates  with  the 
laminas  only  on  one  side  and  over  three  times  those  for  the  aluminum 
plate.  However,  for  increasing  valioes  of  0,  the  values  of  A  decrease 
until  little  or  no  improvement  in  stability  is  obtained  by  using  the 
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conposite  material.  For  values  of  0  between  0°  and  90°,  increasing 
the  number  of  layers  of  composite  material  also  increases  the  flutter 
stability. 

These  results  show  that  composite  materials  may  be  used  very  effec¬ 
tively  to  increase  the  flutter  stability  of  an  isotropic  plate.  For 
square  plates,  aligning  the  fibers  with  the  plate  axis  and  applying  the 
composite  material  symmetrically  about  the  midplane  of  the  plate  both 
result  in  improved  flutter  stability.  If  the  fibers  sure  rotated  with 
respect  to  the  plate  axis,  increasing  the  number  of  laminas  also  in- 
cresuaes  the  flutter  stability. 


Chapter  5 

CONCLUDIHG  REMARKS  AND  RECOMMEHDATIOHS 

For  general  laminated  plates,  the  bending  and  extensional  governing 
equations  are  coupled  and  have  cross-stiffness  terms  which  are  not  in¬ 
cluded  in  classical  orthotropic  plate  theory.  The  coupling  euid  cross¬ 
stiffness  terms  occur  as  a  result  of  the  lamina  principeLl  directions  not 
coinciding  with  the  plate  axis.  These  additional  terms  and  the  coupling 
increase  the  diffictilty  of  obtaining  solutions.  However,  a  solution 
proced\ire  has  been  developed  using  linear  small  deflection  theory  for 
the  flutter  of  arbitrarily  laminated  simply  supported  plates.  The  ex¬ 
tended  Galerkin  method  is  used  to  obtain  solutions  to  the  governing 
equations,  euid  the  aerodynamic  pressure  loading  used  in  the  analysis  is 
that  given  by  linear  piston  theory  with  flow  at  arbitrary  cross-flow 
angles. 

Flutter  solutions  were  obtained  for  typical  symmetric,  angle-ply, 
and  general  laminated  composite  plates,  cud  a  limited  pcurametric  study 
was  conducted.  The  parameters  studied  include  the  number,  orientation, 
and  orthotropy  of  the  lamina;  the  plate  length-width  ratio;  the  inplane 
normal  cud  shear  loads;  cud  the  cross-flow  angle.  In  addition,  flutter 
solutions  for  several  composite  stiffened  aluminum  plate  designs  were 
obtained  to  determine  the  most  flutter  resistant  design. 

The  bending-extensional  coupling  cud  the  cross-stiffhess  terms  both 
have  a  large  destabilizing  effect  on  the  flutter  of  unstressed  laminated 
plates,  but  Increasing  the  number  of  lamincu,  reducing  the  lamina 
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orthotropy,  and  stacking  the  laminas  in  the  "best"  order  reduce  the  de¬ 
stabilizing  effect.  For  a  square  plate,  aligning  the  fibers  with  the 
x-axis  results  in  the  highest  flutter  stability,  but  for  a  plate  with  a 
length-width  ratio  of  2.0,  large  improvements  in  flutter  stability  may 
be  obtained  by  rotating  the  fibers  away  from  the  x-axis.  For  angle-ply 
plates,  inpleuie  normal  and  shear  loads  and  cross- flow  have  a  destabiliz¬ 
ing  effect  on  flutter  similar  to  that  obtained  for  orthotropic  plates. 
However,  for  symmetric  plates  with  the  fibers  not  aligned  with  the 
x-axis,  the  cross-stiffness  terms  result  in  small  inplane  shear  loads 
and  cross-flow  angles  improving  the  flutter  stability.  Flutter  calcula¬ 
tions  for  equivalent  symmetric,  angle-ply,  and  general  unsymmetric 
plates  indicate  that  for  no  cross- flow  and  no  inplane  shear  loads, 
plates  with  an  angle-ply  construction  will  have  the  highest  flutter 
stability.  If  cross-flow  or  inplane  shear  loads  are  present,  syimnetric- 
ally  constructed  plates  may  have  higher  flutter  stability. 

Since  classical  plate  theory  does  not  consider  bending-extensional 
coupling  and  cross-stiffness  terms,  it  gives  inaccurate  and  usually  non¬ 
conservative  flutter  boundaries  for  laminated  plates.  Reduced  bending 
stiffness  theory,  an  approximate  flutter  theory  which  accounts  for  the 
coupling  by  reducing  the  plate  bending  stiffness  as  determined  by  the 
coupling  terms  and  then  neglects  the  coupling  in  solving  the  equations, 
gives  flutter  solutions  that  are  adequate  for  all  plates  for  which  nvimer- 
ical  resxilts  were  obtained.  For  angle-ply  plates,  reduced  bending  stiff¬ 
ness  theory  results  are  obtained  using  published  classical  orthotropic 
plate  theory  solutions.  However,  for  symmetric  or  general  laminated 
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plates,  a  solution  procedure  similar  to  the  present  one  for  symmetric 
plates  is  necessary  to  obtain  even  approximate  solutions. 

The  flutter  stability  of  composite  stiffened  aluminum  plates  was 
considerably  better  than  the  flutter  stability  for  an  equal  mass  aluminum 
plate.  Applying  the  composite  material  to  both  sides  of  the  aluminum 
plate  results  in  better  flutter  stability  than  by  applying  all  the 
material  to  one  side  of  the  plate.  Also,  if  the  fiber  directions  do  not 
coincide  with  the  plate  axis,  increasing  the  number  of  layers  of  material 
improves  the  flutter  stability. 

Since  only  a  limited  parametric  study  was  conducted  in  this  investi¬ 
gation,  it  wovild  be  beneficial  to  use  the  present  analysis  to  conduct  a 
parametric  study  in  greater  depth.  Specifically,  a  wider  variety  of 
lamina  material  properties  should  be  considered  and  additional  plates 
with  general  unsymmetric  construction  should  be  investigated.  Also,  the 
possible  benefits  of  using  composite  materials  to  stiffen  conventional 
orthotropic  and  isotropic  plates  should  be  studied  for  several  additional 
plates.  For  all  results  presented,  the  fiber  directions  of  the  laminas 
were  rotated  by  the  same  absolute  angle  about  the  plate  axis.  Thus, 
further  investigations  should  be  made  into  the  flutter  characteristics 
for  laminated  plates  in  which  the  fiber  directions  of  the  individual 
lamina  are  rotated  independently.  This  becomes  of  special  significance 
for  laminated  plates  with  Inplane  shear  or  cross-flow  where  significant 
iii5)rovement8  in  flutter  stability  may  be  obtained  by  rotating  the  fibers 
away  from  the  plate  axis. 

The  flutter  resvilts  presented  in  the  present  analysis  were  obtained 
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for  simply  supported  plates  with  free  normal  and  tangential  inplane  dis¬ 
placements.  Other  boundary  conditions,  which  restrict  the  inplane  dis¬ 
placements  have  been  considered  in  the  literature  for  simply  supported 
plates  and  have  resulted  in  large  changes  in  natural  frequencies.  Since 
changes  in  flutter  characteristics  are  usually  Msociated  with  changes 
in  the  natureil  frequencies,  the  inplane  boundary  conditions  may  have  a 
significant  effect  on  the  flutter  characteristics.  Thus,  an  extension 
of  the  present  analysis  to  consider  other  inplane  boundary  conditions 
would  be  a  worthwhile  endeavor. 

Although  th-'.  reduced  bending  stiffness  theory  adequately  predicted 
the  flutter  characteristics  for  all  the  plates  studied,  it  was  shown  in 
the  literature  not  to  give  natural  frequencies  and  deflections  as 
accurately  as  desired  for  certain  plates.  Thus,  additional  comparisons 
between  reduced  bending  stiffness  theory  and  the  present  analysis  are 
needed  for  a  wider  range  of  plates.  Also,  since  the  reduced  bending 
stiffness  theory  does  not  account  for  the  inplane  boundary  conditions, 
any  cheinge  in  the  flutter  cheuracteristics  due  to  the  changes  in  boundary 
conditions  suggested  in  the  above  paragraph  would  not  be  shown  with  the 
reduced  bending  stiffness  theory.  These  additional  studies  are  needed 
before  the  reduced  bending  stiffness  theory  can  be  used  with  confidence. 

Since  the  literature  survey  revealed  no  numerical  results  available 
for  general  laminated  plates,  a  significant  contribution  to  the  litera¬ 
ture  could  be  made  by  using  the  present  analysis  to  calculate  the  natural 
vibration  frequencies  and  the  inplane  normal  and  shear  static  buckling 
loads  for  general  laminated  plates.  Althou^  this  study  has  been 
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concerned  primarily  with  flutter  characteristics,  the  analysis  and  re¬ 
sulting  con5)uter  program  can  be  used  without  modifications  to  perform 
the  indicated  calculations. 
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TABLE  I 

LAMINA  MATERIAL  PROPERTHS  USED  IN  ANALYSIS 


Material 

Properties 

Boron-epoxy 

Glass-epoxy 

I-* 

7,500,000  psi 

CM 

CM 

3*000,000  psi 

,‘2,500,000  psi 

®12 

1,000,000  psi 

1,000,000  psi 

CM 

.30 

.30 
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TABLE  II 


MATERIAL  PROPERTIES  FOR  COMPOSITE  STIFFBBED 
ALUMINUM  PLATE 


Material 

Properties 

Aluminum 

Composite 

=11 

29,600,000  psi 

=22 

It .130, 000  psi 

CM 

3«650,000  psi 

2,020,000  psi 

'^12 

.30 

.21 
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TABLE  III 


MATERIAL  PROPERTIES  OF  REFERENCE  PLATE 


Material 

Properties 

Outside 

layers 

Middle 

layers 

^11 

E 

lOE 

CVl 

CM 

lOE 

E 

^12 

.5E 

.5E 

''l2 

.O3U9 

.3^+9 
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TABLE  IV 


COMPARISON  OF  POUR  LOWEST  NATURAL  PRF''rENCIES 


FOR  A  SQUARE  ANGLE-PLY  PLATE  ®12^®22  “  '*12  “ 


Y/E22^^ 

6 

Extended  Galerkin  Method 

0 

19.232 

24.524 

36.229 

54.691 

10 

16.615 

23.584 

36.471 

56.045 

20 

14.610 

24.567 

40.984 

46.414 

30 

14.377 

28.170 

40.452 

49.336 

UO 

14.592 

32.275 

36.235 

58.369 

Rayleigh-Ritz  Method,  ref.  15 

0 

19.232 

24.524 

36.229 

54.691 

10 

16.615 

23.584 

36.471 

56.045 

20 

14.610 

24,567 

40.984 

46.4l4 

^0 

14.377 

28.170 

40.452 

49.336 

UO 

14.592 

32.278 

36.235 

58.369 
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TABLE  V 


COMPARISON  OF  FOUR  LOWEST  NATORAL  FREQUENCIES 
FOR  A  SQUARE  CROSS-PLY  PLATE  »  -5,  ^ 


^11^^22 

Extended  Galerkin  Method 

1 

6.155 

15.148 

15.148 

24 . 619 

10 

7.932 

21.274 

21.274 

31.727 

20 

9.152 

25.136 

25.136 

36.608 

30 

10.210 

28.403 

28.403 

4o.84o 

UO 

11.164 

31.312 

31.312 

44.656 

Rayleigb-Rltz  Method, 

ref.  15 

1 

6.155 

mmm 

15.148 

24.619 

10 

7.931 

21.275 

21.275 

31.721 

20 

9.152 

25.135 

25.135 

36.601 

30 

10.210 

28.400 

28.400 

40.831 

UO 

11.163 

31.308 

31.308 

44.646 
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Figure  1.-  Geometry  of  angle-ply  plates. 


Figure  2.-  Geometry  of  cross-ply  plates. 
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Figure  3.--  Lamina  coordinate  systea. 
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z,3 


Figure  Rotation  of  lamina  axes 
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Figure  6.-  Positive  direction  of  stress  and  moment  resultsuits. 
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90 


Plate  P-l 


Figure  10.-  Nomenclature  and  geometry  of  general  laminated  plates 
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Figure  12.-  Convergence  of  solution  for  a  symmetric,  boron-epoxy, 
plate  with  K  »  U. 
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Cb)  0  =  60°. 


Figure  13.-  Concluded. 


Figure  lU.-  Comparison  of  flutter  solutions  for  square  reference  plate 
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Figure  15.“  Converged  flutter  boimdaries  for  sqTiare  reference  plate. 
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101 


Cb)  Boron-epoxy  plate. 
Figure  l6.-  Concluded. 
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Figure  l8.-  Flutter  boundaries^  for  squ^e,  boron-epoa^,  Bynnetric  plate 
with  Inplane  normal  loads.  n  ■*9  •>  *  4, 

y 


Figure  19. «  Flutter  boundaries  for  square,  boron-epoxy,  symmetric  plate 
with  inplane  shear  loads,  ^ 


9,  deg 


Figure  23.-  Flutter  boundaries  for  boron-epoxy,  auagle-ply  plate  with 
a/b  =  2.0.  '  0* 


A.  deg 


Figure  25.-  Flutter  boundaries  for  sq.uare,  boron-epoxy,  angle-ply  plate 
with  cross-flow.  =*  *  N _ =  0;  K  =  1+. 
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(a)  Plate  P-1. 


Figure  26.-  Flutter  boundeuries  for  square,  boron-epoxy,  general  laminated 


plates.  N  =  N 
X  y 


N  =  0. 
xy 


Figure  27.-  Compeurison  of  flutter  boundaries  for  square,  symmetric 

angle-ply,  and  genersa  laminated  plates.  N„  =  N  = 

y 


